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ABSTRACT 

We present here the first study of the X-ray properties of an evolutionary sample 
of merging galaxies. Both ROSAT PSPC and HRI data are presented for a sample of 
eight interacting galaxy systems, each believed to involve a similar encounter between 
two spiral discs of approximately equal size. The mergers span a large range in age, 
from completely detached to fully merged systems. 

A great deal of interesting X-ray structure is seen, and the X-ray properties of each 
individual system are discussed in detail. Along the merging sequence, several trends 
are evident: in the case of several of the infrared bright systems, the diffuse emission 
is very extended, and appears to arise from material ejected from the galaxies. The 
onset of this process seems to occur very soon after the galaxies first encounter one 
another, and these ejections soon evolve into distorted flows. More massive extensions 
(perhaps involving up to 10 10 M & of hot gas) are seen at the 'ultraluminous' peak of 
the interaction, as the galactic nuclei coalesce. 

The amplitude of the evolution of the X-ray emission through a merger is markedly 
different from that of the infrared and radio emission however. Although the X-ray 
luminosity rises and falls along the sequence, the factor by which the X-ray luminosity 
increases, relative to the optical, appears to be only about a tenth of that seen in the 
far-infrared. This we believe, may well be linked with the large extensions of hot gas 
observed. 

The late, relaxed remnants, appear relatively devoid of gas, and possess an X-ray 
halo very different from that of typical ellipticals, a problem for the 'merger hypothe- 
sis', whereby the merger of two disc galaxies results in an elliptical galaxy. However, 
these systems are still relatively young in terms of total merger lifetime, and they may 
still have a few Gyr of evolution to go through, before they resemble typical elliptical 
galaxies. 

Key words: galaxies: individual: (Arp 270, Arp 242, NGC 4038/9, NGC 520, Arp 220, 
NGC 2623, NGC 7252, AM 1146-270) - galaxies: interactions - galaxies: evolution - 
galaxies: ISM - galaxies: peculiar - X-rays: galaxies 



Galaxies were once thought of as 'island universes', evolving 
slowly in complete isolation. This is now known not to be 
the case. Galaxies interact in a variety of ways with their 
environment, both with satellite and neighbouring galax- 
ies, and, in the case of galaxies within groups and clusters, 
with large masses of tenuous hot gas. Collisions and merg- 
ers of galaxies are now thought to be one of the most dom- 
inant evolutionary mechanisms (Schweizer 1989). Indeed, if 
events during galaxy formation are counted, there are prob- 
ably very few galaxies that were not shaped by interactions 
or even outright mergers (Toomre 1977). The position of a 
galaxy in Hubble's (1936) morphological sequence may in 



1 INTRODUCTION 



fact depend mainly on the number and severity of merger 
events in its past history. Pure disc systems, formed from 
relatively isolated protogalactic gas clouds, appear at one 
end of the Hubble sequence, the giant ellipticals, possibly 
produced through mergers of similar spirals, appear at the 
other, and in between, mergers between galaxies of differing 
mass produce galaxies with a wide range of bulge to disc 
ratios. 



This merger hypothesis discussed above, the idea that 
elliptical galaxies might be formed from the merger of two 
disc galaxies, was first suggested by Toomre (1977). The 
conversion of orbital to internal energy during a close tidal 
encounter causes the two progenitor systems to sink together 
and coalesce violently into a centrally condensed system, dis- 
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rupting any pre-existing discs, and largely randomizing the 
stellar motions, as shown in many N-body simulations (e.g. 
Barnes 1988). The problem that ellipticals have central mass 
densities some 10 2 — 10 3 times higher than normal spirals 
has been overcome by the discovery that young merger rem- 
nants contain huge central concentrations of molecular gas 
(Sanders et al. 1986; Sanders et al. 1988a; Scoville & Soifer 
1991). Also, elliptical-like luminosity profiles have been dis- 
covered in the infra-red underlying the disturbed optical 
morphology of merging galaxies (Wright et al. 1990). There 
have been many attempts to disprove the merger hypoth- 
esis. Counter-evidence cited includes the fractional abun- 
dances of ellipticals in clusters, the colour-luminosity and 
metallicity-luminosity relationships, radial gradients within 
ellipticals and with the incidence of globular clusters within 
ellipticals (see Hernquist 1993). These criticisms can mostly 
be overcome by appealing to a variety of ideas, such as the 
growth of substructure in hierarchical cosmogonies, the in- 
complete nature of violent relaxation, and merger induced 
star-formation. 

Studies of the evolution of merging galaxies have been 
performed in many regions of the electromagnetic spectrum 
(Joseph & Wright 1985; Telesco, Wolstencroft & Done 1988; 
Casoli et al. 1991; Hibbard & van Gorkom 1996), and these 
are discussed in detail in both Section |^ and Section ^. One 
region of the electromagnetic spectrum that has been ne- 
glected so far though, has been the X-ray. 

The ROSAT X-ray telescope (XRT), with the Position 
Sensitive Proportional Counter (PSPC) (Pfeffermann et al. 
1986) at the focal plane, offers three important improve- 
ments over previous X-ray imaging instruments, such as the 
Einstein IPC. First of all, its spatial resolution is very much 
improved, the 90% enclosed energy radius at 1 keV being 
27" (Hasinger et al. 1992). Secondly, the spectral resolution 
of the PSPC is significantly better (AE/E ~ 0.4 FWHM at 
1 keV) than earlier X-ray imaging instruments and this al- 
lows the derivation of characteristic source and diffuse emis- 
sion temperatures. Finally, the PSPC internal background is 
very low (~ 3 x 10~ 5 ct s" 1 arcmin ; Snowden et al. 1994) , 
thus allowing the mapping of low surface brightness emis- 
sion. The High Resolution Imager (HRI) on the other hand, 
because of its excellent spatial resolution (~ 5") and relative 
insensitivity to diffuse emission, is an ideal instrument for 
further investigation into the point source populations. 

Furthermore, as ROSATs energy band is relatively soft 
compared to that of the Einstein IPC (and to other previ- 
ous instruments), the hot halos now known to exist around 
many nearby normal galaxies (e.g. Bregman & Pildis 1994; 
Wang et al. 1995; Read, Ponman & Strickland 1997), may 
be observable within these more distant, though more ac- 
tive systems. Indeed as shall be discussed later, a hot halo is 
known to exist around the merging galaxy pair, the Anten- 
nae (NGC 4038/9) (Read, Ponman & Wolstencroft 1995). 

In the ROSAT band, it appears that the X-ray emis- 
sion from normal spiral galaxies is made up of a complex 
mix of sources. Both stellar sources such as low- and high- 
mass X-ray binaries, cataclysmic variables and normal main 
sequence stars, and interstellar sources such as supernova 
remnants (SNRs) and the hot phases of the ISM, contribute 
to the total X-ray emission. The merger of two similar spiral 
galaxies may lead to additional sources of X-ray emission. 
The advances in instrumentation discussed above should al- 



low a substantial improvement in our knowledge of the com- 
plex X-ray properties of interacting galaxies over what was 
possible with Einstein (see Fabbiano 1989 and references 
therein). The variation in X-ray emission with merger stage 
and position can be compared with corresponding multi- 
wavelength results to reveal new clues as to how the tidal 
forces trigger star formation, and to constrain the nature 
and duration of the activity which ensues. 

Studies of the nearest, and therefore best resolved, in- 
teracting systems are especially important, and many of the 
most attractive targets have been observed by ROSAT. We 
have undertaken a programme of ROSAT PSPC and HRI 
observations of a carefully chosen chronological sequence of 
merging galaxies in order to study the evolution of their X- 
ray properties through the merging process. We have com- 
bined our own ROSAT data with other archival ROSAT 
data and with individual studies published recently. Where 
possible, we have attempted to analyse these systems in a 
uniform manner. In particular we have carefully separated 
point sources from what appears to be diffuse X-ray emis- 
sion, and have determined the spectral properties of each. 

In the present paper we describe the analysis we have 
performed, and present the basic results for each of our sam- 
ple systems in turn. In some cases, we are able to compare 
with previous authors' results from the same data, but in 
many instances our analysis is the first reported ROSAT 
analysis, and in some cases it represents the first picture of 
the X-ray emission from an imaging telescope. 

The plan of this paper is as follows. After some com- 
ments (Section |^) on the selection of the sample, Section ^ 
describes the observations and the data reduction methods. 
The results for each galaxy are presented in Section El to- 
gether with notes on the individual systems. In SectionlBL we 
comment on the range of galaxy properties which emerges, 
and discuss the X-ray evolution of the sample as a whole. 
Finally, in Section o we present our conclusions. 



2 THE MERGING GALAXY SAMPLE 

The basic problem with a chronological study is to establish 
the ordering of the sample. The most famous chronological 
sequence of merging galaxies is undoubtedly the 'Toomre se- 
quence' (Toomre 1977), containing the 11 best examples of 
ongoing mergers from the New General Catalogue of Neb- 
ulae and Clusters of Stars (NGC). These systems appear, 
when ordered chronologically, to range from well separated, 
strongly interacting systems, apparently destined to even- 
tually merge, through systems in the throes of merging, to 
well merged systems. This sequence has formed the basis 
for other evolutionary studies made at other wavelengths, 
as discussed below, and indeed, forms the basis for the se- 
lection of the present sample. 

Relevant studies at other wavelengths were, until very 
recently, rare. Joseph and Wright (1985) studied highly dis- 
turbed systems in which the two disc galaxies had lost their 
individual identities and appear as a single coalesced object, 
the primary morphological indication of which is the pres- 
ence of the two tidal tails. Dynamical simulations (Toomre 
and Toomre 1972) suggested that the tails will persist for 
~ 10 9 years, and the faintness of these tails, together with 
the degree of coalescence, were used as indicators of rel- 
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ative merger age. Comparing the infrared luminosities of 
these galaxies they found that the two youngest and two old- 
est systems were less luminous than the four 'middle-aged' 
systems, suggesting that merging galaxies undergo a phase 
of starburst activity of characteristic luminosity ~ 10 12 Lq. 
Another useful indicator of the progress of a starburst was 
found by Telesco, Wolstencroft & Done (1988); for compa- 
rably sized interacting galaxies, the dust colour temperature 
and the interaction strength is greatest for the pairs with the 
smallest separation. Furthermore, Casoli et al. (1991), while 
studying the molecular gas content within a series of merg- 
ing systems, discovered that the far-infrared colour temper- 
ature, the Lfir/Lb ratio and the Lfir/M{Y{2) ratio all 
evolve along the sequence, increasing in the early stages of 
the collision, reaching a climax, and then decreasing to val- 
ues typical of ellipticals. 

Much more recently, Hibbard & van Gorkom (1996) 
have presented Hi, Hn, and R-band observations of an evo- 
lutionary sequence of merging systems, again taken primar- 
ily from the Toomre sequence. Their five systems, four of 
which appear in the Toomre sequence, and three of which 
appear in this present paper, span a large range in merger 
stage, unlike the work of Joseph & Wright (1985) discussed 
above, where only middle stage mergers were considered. 
Much of Hibbard & van Gorkom's (1996) work is very rel- 
evant to the present paper, and many of their results are 
used in our discussion (Section ^). They find that, moving 
from early to late-stage mergers, a larger fraction of the 
neutral hydrogen is found outside of the optical confines of 
the system, and in the last stage, all of the Hi is found 
in the tails, with none in the central regions. They suggest 
that the interaction-induced relaxation causes most of the 
atomic gas within the original discs to compress, condense, 
form stars and/or shock-heat to X-ray temperatures, any 
residual atomic gas remaining in the remnant being heated 
via the resultant starburst. 

In order to make the establishment of a reliable chrono- 
logical ordering as easy as possible, the systems were care- 
fully chosen on the basis of both close attention to the 
work of the previous authors' described above, especially 
Toomre's (1977) sequence, and the following important cri- 
teria: 

• All the systems contain, or appear to have evolved from, 
two spiral galaxies of fairly equal mass. 

• A wide range of merger stages, wider in fact than the 
Toomre sequence, from completely detached to fully merged 
systems, is covered. 

• All the systems are infrared and radio bright, indicating 
the presence of unusual star-formation activity. 

• Low absorbing columns (1.2 — 5.2 x 10 20 cm -2 ) exist in 
the directions of all these systems, maximizing the sensitiv- 
ity to soft X-ray emission. 

• All the systems are large enough for worthwhile spatial 
resolution to be achieved with the PSPC, and for detailed 
mapping of point sources with the HRI. 

• Much additional multiwavelength information is avail- 
able for all these systems. 

Optical morphology (especially the appearance of tidal 
tails, and the proximity of the two nuclei) has been used 
as the primary measure of merger stage, but this has been 
supplemented with the following additional indicators. 



Toomre & Toomre (1972) have shown that proper tail 
construction requires that the perturbing mass be at least 
comparable to the perturbed, and that therefore, major tidal 
tails should protrude from both of two interacting discs of 
roughly equal dimensions as in NGC 4676 (the Mice) and 
NGC 4038/9 (the Antennae). For systems where the central 
parts of the galaxy have not yet merged, tail length seems 
the best indicator of age, increasing with time (Toomre & 
Toomre 1972), although tail length also increases with the 
closeness of the encounter. Once the nuclei of the galaxies 
have merged, the tails begin to diffuse away, and so using 
the morphology of the systems to order them chronologically 
comes down to a combination of tail length, tail faintness 
and degree of coalescence of the nuclei. 

Listed below, in chronological order, are the eight sys- 
tems within the sample, selected on the basis of the above 
considerations. Optical pictures of the eight systems are 
shown in figure |l| All are taken from the digitised sky sur- 
vey, the southern systems from the U.K. Schmidt plates and 
the northern systems from the Palomar Schmidt plates. 

Arp 270 (NGC 3395/6) consists of two low IR-flux 
galaxies, obviously involved in some sort of interaction, 
though the lack of tidal tails is suggestive of them not hav- 
ing properly encountered one another yet. This system is in 
fact a pre- Toomre sequence system. 

Arp 242, the Mice (NGC 4676), occurs second in the 
Toomre sequence and, though it appears rather similar to 
Arp 270, the presence of tidal tails and an obvious bridge 
between the two galaxies, together with an increase in IR 
activity, indicates that the galaxies have begun to interact, 
and have passed each other. 

NGC 4038/9, the Antennae (Arp 244), is a classic 
example of an interacting system with two equally long tails 
and distinct central masses. Although it occurs first in the 
Toomre sequence, CO (Stanford et al. 1990), radio (Hummel 
& van der Hulst 1986) and X-ray observations (Read et al. 
1995) indicate that the discs have begun to merge, and the 
Antennae therefore should really be placed sometime after 
the Mice. 

NGC 520 (Arp 157), seventh in the Toomre sequence, 
and classified as an intermediate-stage merger by Hibbard 
& van Gorkom (1996), is as radio and infrared bright as the 
Antennae, and has two smaller tails as well as two nuclei 
and two velocity systems in its spectra, indicative of a young 
merger. 

Arp220, one of the superluminous IRAS galaxies, is a 
prototypical merger. It is the brightest of the sample across 
the whole of the electromagnetic spectrum (except for the 
optical). Huge Ha filaments exist (Heckman, Armus & Mi- 
ley 1990), suggestive of vigorous starburst activity. Two dis- 
tinct (though very close) nuclei are visible in the infrared 
(Majewski et al. 1993). 

NGC 2623 (Arp 243), appearing eighth in the Toomre 
sequence, is, like Arp 220, a superluminous IRAS galaxy, and 
is also very bright in the radio. Very long tails are visible 
but the central masses have become indistinguishable. Only 
one true nucleus is thought to exist. 

NGC 7252 (Arp 226), the atoms for peace galaxy, 
is the prototypical merger remnant, and is the last sys- 
tem within the Toomre sequence. It has two equally long, 
very faint tails emanating from a disturbed, double motion 
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spheroidal body whose mean light distribution closely fol- 
lows the de Vaucouleurs r 1//4 law, typical of ellipticals. 

AM 1146-270, beyond the end of the Toomre se- 
quence, is a spheroidal galaxy with a small tail-like struc- 
ture, appears to be the site of a great deal of recent star 
formation, hence the many blue knots surrounding its nu- 
cleus. 

Table jj] lists the sample systems together with the fol- 
lowing basic properties. Distances are, where possible, taken 
from Tully (1988). These are based on H = 75 km s" 1 
Mpc - , and assume that the Galaxy is retarded by 300 km 
s" 1 from universal expansion by the mass of the Virgo clus- 
ter. For the more distant galaxies, distance values are ei- 
ther taken from Condon et al. (1990), where the same Ho 
and Virgocentric flow correction is used, or are calculated 
from the group-flow-corrected radial velocities given in the 
Second Reference Catalog of Bright Galaxies (RCBG) (de 
Vaucouleurs, de Vaucouleurs & Corwin 1976), using a Ho 
of 75 km s -1 Mpc" 1 . AM 1146-270's radial velocity is taken 
from Sekiguchi & Wolstencroft (1993). 

Optical (B) luminosities for all the systems in the sam- 
ple are calculated as in Tully (1988), i.e. as a re-expression 
of the absolute magnitude, itself following from the blue ap- 
parent magnitude Bt, and the distance D: 

logL fl (L ) = 12.192 — OABt + 2 log D. 

Blue apparent magnitudes are, again where possible, 
taken from Tully (1988). Values for the more distant sys- 
tems are taken from RCBG (de Vaucouleurs, de Vaucouleurs 
& Corwin 1976) and from NGC 2000.0 (Dreyer 1988) (the 
value for AM 1146-270 is taken from Smith & Hintzen 
1991)). 

FIR luminosities are calculated from IRAS 60 and 
100 \im fluxes (taken from the IRAS Point Source Cata- 
logue) using the expression 

Lfir — 3.65 x 10 5 [2.58S , 60nm + Swo^m] D 2 Lq, 

(e.g. Devereux & Eales 1989). Here D is the distance 
in Mpc and Seo/j.m. and Sioo/am are the IRAS 60 and 100 fim 
fluxes (in Janskys) . Also given in table ^ are the infrared to 
blue luminosity ratio, Lfir/ Lb, and the dust colour tem- 
perature, Seo/Sioo- These two infrared indicators are the 
two suggested by previous studies (Joseph & Wright 1985; 
Telesco et al. 1988) to be worthwhile indicators of starburst 
activity, and it can be seen that they rise and fall along the 
chronological sequence in just the way one might expect if 
a burst of star formation were trig gered during a merger. 
Finally, included also in Table [12| are the radio luminosi- 
ties from Condon et al.'s (1990) 1.49 GHz atlas of the IRAS 
bright galaxy sample, supplemented in the cases of Arp 242 
and NGC 7252 (upper limit), with the luminosities (con- 
verted to Ho = 75km s -1 Mpc -1 ) from Heckman's (1983) 
study. 



3 OBSERVATIONS AND DATA REDUCTION 

Both PSPC and HRI data have been analysed within this 
study, and the methods of data extraction and reduction for 
both types of data are described below. All of the HRI data 
and the PSPC data for three of the systems (NGC 4038/9, 
NGC 2623 and AM 1146-270) were obtained through our 



own pointed observations. The remaining PSPC datasets 
were obtained from the UK ROSAT Data Archive Centre at 
the Department of Physics and Astronomy, Leicester Uni- 
versity U.K. 

3.1 PSPC Observations and Data Reduction 

All of the PSPC datasets have been analysed in essentially 
the same way using the STARLINK ASTERIX X-ray anal- 
ysis system. This method is described below, and in greater 
detail in Read et al. (1997). Departures from this standard 
procedure are described in detail in the notes on the indi- 
vidual systems (sections [O— |l.8|). In every case, the target 
system in question lay at the centre of the PSPC field, thus 
minimizing the effects of vignetting and blurring. 

Once the data were 'cleaned' of high background pe- 
riods (« 2 — 3% of the data), they were binned into a 
(0.1— 2.3keV) spectral image (or 'data cube') of size ap- 
proximately twice that of each system. An annulus situated 
outside the PSPC central support ring, and with the bright 
sources within it removed, was formed from the data, and 
this background region was, together with knowledge of the 
PSPC vignetting function, used to construct a background 
model. When this was subtracted from the data cube, the 
resultant background-subtracted data cube could be used to 
form images in different spectral bands. 

PSS (Allan, Ponman & Jeffries, in preparation), a point 
source search program which uses a likelihood technique to 
search for enhancements above the background, was used 
to search for point-like emission within these images. The 
basic method of PSS involves a comparison of the input 
dataset (in this case, the image) with a model, comprising 
of a background and a scaled PSF. The positions of the 
detected point sources, i.e. those with a significance > 4<r 
(except where stated in the text), were cross-correlated with 
a variety of stellar and non-stellar catalogues, including the 
SIMBAD catalogue and the Einstein and IRAS point source 
catalogues. 

As in Read et al. (1997), we attempted to separate the 
diffuse from the source emission by removing data from the 
background-subtracted data cube at the position of each 
PSS source over a circle enclosing 93% of the energy of a 
0.5 keV point source (as all the sample systems subtended 
small angles, this extraction radius was never greater than 
about half an arcminute) . The remaining data were then col- 
lapsed into a spectrum and corrected for vignetting effects 
and exposure time. To account for the diffuse flux lost in 
the source removal process, the diffuse spectrum was renor- 
malised using a 'patched' image, where the 'holes' left after 
source removal were filled by bilinear interpolation. 

This method of separating the point source emission 
from the diffuse emission did not, in general, work as well as 
in Read et al. (1997). This is not surprising as the systems in 
the present study are at much greater distances than those 
studied in Read et al. (1997), and both limitations on sen- 
sitivity and on resolution have had a more serious impact. 

As in Read et al. (1997), source spectra were binned 
directly from the raw data from circles of the same size as 
those used in the source removal procedure above. To ob- 
tain true source spectra however, both the local true back- 
ground flux, obtained by scaling the background emission 
from source-free regions of the field using knowledge of the 



© 0000 RAS, MNRAS 000, 



The X-ray Evolution of Merging Galaxies 5 



Figure 1. The merging galaxy sample: A) Arp 270, B) Arp 242, C) NGC 4038/9, D) NGC 520, E) Arp 220, F) NGC 2623, G) NGC 7252, 
H) AM 1146-270 



System Other names Distance log Lb log LpjR Lpjr/Lb ■S6o/>S i ioo logL rat j 

(Mpc) (ergs" 1 ) (ergs -1 ) Whz -1 



Arp 270 


NGC 3395/6 


28 


44.08 


43.61 


0.335 


0.486 


21.93 


Arp 242 


NGC 4676 


88 


44.06 


44.09 


1.069 


0.534 


22.18 


NGC 4038/9 


Arp 244 


25 


44.19 


44.20 


1.026 


0.516 


22.64 


NGC 520 


Arp 157 


28 


43.87 


44.15 


1.900 


0.651 


22.24 


Arp 220 


UGC 9913 


76 


43.96 


45.52 


36.04 


0.884 


23.39 


NGC 2623 


Arp 243 


78 


43.78 


44.87 


12.28 


0.862 


22.85 


NGC 7252 


Arp 226 


63 


44.28 


44.00 


0.526 


0.565 


<22.78 


AM 1146-270 




25 


42.93 


42.55 


0.416 


0.590 





Table 1. The merging galaxy sample, listed in order of evolutionary age. Sources for the tabulated parameters are given in the text. 



PSPC vignetting function, and the local diffuse flux, cal- 
culated by scaling the integrated spectrum of the diffuse 
emission with an estimate of the fraction of the diffuse flux 
underlying each source, had to be removed from each source 
spectrum. 

So, in summary, three types of spectrum were formed 
from each dataset: 



• An integrated Spectrum was extracted from a cir- 
cular area of diameter slightly greater than the optical di- 
ameter of each system. Although particle and cosmic X- 
ray backgrounds were removed, no sources were excluded. 
Standard spectral models (power law, bremsstrahlung and 
Raymond and Smith (1977) hot plasma) were used in the 
fitting. Even though these spectra contain emission from 
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everything around each system's position, including galac- 
tic sources, galactic diffuse emission, and foreground and 
background sources (both Galactic and extragalactic), some 
of the resultant spectral fits are very acceptable. The net 
number of counts, the \ 2 values, number of degrees of 
freedom (n.d.o.f.), and the integrated X-ray luminosities 
(0.1— 2.0keV) inferred from the spectral fitting, are shown 
in each of Tables | to [to| 

• A Diffuse spectrum was extracted, using the method 
described above, from an area the same size as that 
used for the integrated spectrum. The correctly normal- 
ized (0.1 — 2.0keV) X-ray luminosities of each system's dif- 
fuse emission are tabulated in Tables [] to along with the 
net diffuse counts and the results of the spectral fits (us- 
ing again the same models as those used above) . Because of 
low count rates, x 2 fitting could only be used in the case of 
the Antennae. In all other cases, a maximum likelihood cri- 
terion, which allows for the Poissonian nature of the data, 
had to be used. Unfortunately, when using likelihood fit- 
ting, the absolute value of the statistic is not informative 
(and hence we have not given it). The Cash statistic (Cash 
1979) can be used however, to compare the relative quality 
of different fits, and to derive confidence intervals for model 
parameters, since differences in the value of the Cash statis- 
tic from one model to another, are distributed in the same 
way as % 2 ■ Furthermore, only the freezing of certain param- 
eters at sensible values led to acceptable fits in a few cases 
(these are indicated by a bracket ted 'F'). All these results 
are described in detail in Sections 4.1 to 4. £. Because of the 
large distances involved, many of these systems' diffuse emis- 
sion spectra will be complicated, involving some genuinely 
diffuse gaseous emission, and some contribution from both 
unresolved point sources and from 'stumps' left in the sub- 
traction of the bright sources. Hence, we attempted, where 
possible, to fit a 2-component model to the diffuse spectrum, 
comprising of a Raymond and Smith hot plasma (represent- 
ing the truly diffuse gas) and a highly absorbed hot (10 keV) 
component (representing unresolved sources). These results 
are also described in more detail later in the notes on the 
individual systems. 

• Individual source spectra were extracted as de- 
scribed above. Sources known not to be connected with the 
galaxy were excluded from the analysis, as were sources that 
lay outside both the optical emission and any diffuse X- 
ray emission. The results of the spectral analysis of all the 
sources within each system are tabulated below in each of 
the individual galaxy sections. The point source spectra gen- 
erally contained low numbers of counts, and Gaussian statis- 
tics could not be assumed. A maximum likelihood criterion, 
as described above, was used. 

3.2 HRI Observations and Data Reduction 

The three HRI datasets (NGC 520, NGC 2623 and 
NGC 7252) were analysed in exactly the same way, again 
using the STARLINK ASTERIX X-ray analysis system. In 
each case the data were first binned into a 0.4° x 0.4° im- 
age with a 3" resolution, thus exploiting the HRI's superior 
spatial resolution. Background subtraction of HRI data, al- 
though significantly easier than of PSPC data, is still non- 
trivial, and two methods were utilised, which were compared 
both with each other and with the mean HRI background 



System 


BR 


BR 




0.3 kcV 


3.0 keV 


NGC 520 


6.47xl0~ 4 


3.50X10" 4 


NGC 2623 


4.54X10" 4 


2.31xl0~ 4 


NGC 7252 


7.11X10" 4 


4.99xl0~ 4 



Table 2. Counts-to-flux Conversion factors for the ROSAT HRI 
detector for the three individual HRI datasets and for two sep- 
arate thermal bremsstrahlung models. The values are corrected 
for Galactic absorption and are in units of 10 — 11 erg cm -2 s~ x 
cts" 1 . 



rate (» 1.1 x 10~ 6 ct s _1 arcsec" 2 ; David et al. 1993). Firstly, 
a polar profile of the data was formed about the centre 
of the system, with radial bins of width 0.001°. This was 
seen to level out at large radius to a constant background 
level. The second method involved an evaluation of the mean 
value, once all the bright features had been removed from the 
data. The background values obtained from the two methods 
agreed well in all three cases. These values were subtracted 
from the raw data (over such a small area, the vignetting is 
essentially flat, and can be ignored). 

The resultant image was then exposure corrected, and 
searched for sources, again using PSS. Again, a detection 
threshold of 4a was assumed. Images, smoothed on many 
different spatial scales (from 3"to 18"), were formed to 
search for extended low-surface brightness features. 

Count rates of the detected point sources within the 
three HRI datasets were converted into fluxes (and then 
into luminosities), assuming 0.3 keV and 3.0 keV thermal 
bremsstrahlung models. HRI conversion factors for both 
models from number of counts into fluxes (taking into ac- 
count the different exposure times and Galactic foreground 
Nh values) are given in Table ^. 

Although the spectral response of the HRI is rather 
poor, it is thought (Fraser 1992) that the ratio of the num- 
ber of counts in channels 1—5 to the number of counts in 
channels 6—11, used by Wilson et al. (1992), does seem the 
most sensitive energy indicator that can be constructed. The 
above analysis was also repeated for data within these two 
energy bands. 



4 RESULTS 

The following eight subsections describe the results of the 
analyses of ROSAT PSPC and, where applicable, HRI data 
of the sample systems. Where the analysis technique differs 
from that described in Section ^ this is noted in the text. 

Within each of Sections to |4.8[ a table is given, show- 
ing the main results of the analysis of each set of ROSAT 
PSPC data. The results of the spectral analyses of each of 
the three types of spectra formed as described in Section ^ 
(an integrated spectrum, a diffuse spectrum and individual 
source spectra), are given. The number of counts, plus the 
statistical uncertainty, contained within each spectrum are 
given, together with the results of the best spectral fit to 
each spectrum. A quoted photon index indicates that the 
best fit is a power law model, whereas if instead, a temper- 
ature is given, the best fit is either a Raymond & Smith 
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hot plasma model (if the metallicity is quoted), or it is a 
bremsstrahlung model (if no metallicity is quoted). Where 
X 2 -fitting has been used (in fitting most of the integrated 
spectra plus a small number of other high-count rate cases) , 
the values of the % 2 statistic are given together with the 
number of degrees of freedom. The (0.1— 2.0keV) X-ray lu- 
minosities tabulated are those escaping from the system in 
question, i.e. before absorption in our own galaxy. Also given 
in each case is the value of Lx(diffuse-l-sources), the sum 
of the diffuse X-ray luminosity and each of the individual 
source luminosities. 

Unless otherwise stated, we also present for each source 
a PSPC image showing contours of background-subtracted, 
exposure-corrected X-ray emission in the 0.1—2.3 keV band, 
superimposed on an optical image of the system. These X- 
ray images have a resolution of 5" and have been lightly 
smoothed with a Gaussian of standard deviation 10" to sup- 
press noise. The contour levels increase by factors of two 
from 7.2xl0~ 4 ct s _1 arcmin -2 . 



4.1 Arp 270 

Arp 270 (also VV 246; Vorontsov-Velyaminov, 1959) is com- 
posed of two galaxies, thought to be of comparable masses 
(Davis & Seaquist 1983), NGC 3395, to the west, and the 
somewhat smaller, more irregular galaxy, NGC 3396, to the 
east. The angular separation of their nuclei (« 90") corre- 
sponds to a separation of about 12 kpc at the distance of 
28 Mpc assumed within this paper. In the far-infrared, they 
appear as a strong, though unresolved source (Soifer et al. 
1987), and in the UV, they are noted in the KISO cata- 
logue (Takase & Miyauchi-Isobe 1987) as an interesting pair 
of UV-excess galaxies. Many radio continuum studies have 
been made, the most noteworthy being that of Huang et 
al. (1994), who find that the radio emission in each galaxy 
comes predominantly from that part of each galaxy near its 
companion galaxy. Whereas NGC 3396 shows two strong ra- 
dio peaks and one weak one, its partner has a much more 
complex structure, with six or more peaks in the north- 
eastern part of the galaxy alone. A clear bridge between 
the two galaxies is also seen in the radio, roughly coincident 
with the optical. 

Arp 270 has been observed in X-rays previously to 
ROSAT, with the Einstein IPC. The Einstein image (Fab- 
biano et al. 1992) shows very little structure, what there is 
appearing to be centred on NGC 3395. Suggestions of ex- 
tensions are also visible to the north-west and south-west of 
the system. 

The ROSAT image (Fig. ^) shows far more structure 
than seen with Einstein. Emission is seen to come from both 
galaxies, and this emission appears to be enshrouded in more 
diffuse emission. The X-ray feature associated with the west- 
ern galaxy (source A; see Table ^ appears to be centrally 
positioned and lies less than 2" from the position given in 
the Third Reference Catalog of Bright Galaxies (RC3) (de 
Vaucouleurs et al. 1991). Source B however, does appear, 
like the radio emission discussed above, to emanate from a 
position in NGC 3395 nearer to its companion galaxy (X-ray 
source B lies some 27" from the RC3 position of NGC 3396. 
As Huang et al. (1994) point out, this is not expected from 
dynamical simulations. Instead, the gas should be forced 
through the tidal interaction towards the nuclear regions of 



Figure 2. Contours of (0.1— 2.3keV) X-ray emission shown su- 
perimposed on an optical image of Arp 270. Contour levels in- 
crease by factors of two from 7.2 X 10 — 4 ct s" 1 arcmin~ 2 .) 



each galaxy, and be compressed into star-formation. This 
process however, as Huang et al. note, will take some time, 
and the star-formation may well begin before the gas reaches 
the nuclear regions. Non-nuclear bursts of star-formation 
can occur in the collision between two gas-rich spiral galax- 
ies, when a burst of massive star-formation is first triggered 
in the overlapping/interpenetrating regions of the galaxies 
(Jog & Solomon 1992). This overlapping region will at first 
be the very outer discs of the two galaxies. In the case of 
Arp 270 though, it appears that this very initial phase may 
be nearly over, as gas is now seen close to the nuclear re- 
gions, i.e. the forcing of gas through the interaction towards 
the galactic nuclei has begun. The fact that, in the case of 
NGC 3396, this hot gas is seen only close to the nucleus, 
not yet within it, may indicate that this second phase, the 
tidal forcing of gas towards the two nuclei, may have only 
just begun. Furthermore, there is quite a good correlation 
between the X-ray contours in Fig. ^| and the radio data of 
Huang et al. (1994) - both the tentative X-ray 'bridge', seen 
apparently connecting the two galaxies, and the radio bridge 
appear to connect the two systems along the southern edge 
of the optical bridge. 

Both source spectra are best fit with low-temperature, 
absorbed plasma models, suggestive of the emission from 
both regions being due to hot gas, rather than to evolved 
stellar components. The fact that no evidence was found 
in the radio for an active galactic nucleus in either galaxy, 
and that the radio spectra of the galaxies are relatively flat 
(at low frequencies) (Huang et al. 1994) , indicative of a good 
deal of supernova activity, adds credence to the hot gas idea. 

The two other detected sources in the field, the strong 
point-like source to the north of NGC 3396, and the weaker 
source within the X-ray 'tail' to the south-west of NGC 3396, 
appear to explain the suggested extensions tentatively seen 
in the Einstein image. The northern source is unlikely to 
have anything to do with the system, and is most probably 
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a background quasar, its best spectral fit being a power-law 
fit with a photon index, a = 2.3 ± 0.5, consistent with the 
mean spectrum of quasars in the ROSAT band (2.2±0.2; 
Branduardi-Raymont et al. f994; Roche et al. 1995). The 
X-ray feature to the south-west of NGC 3395 however, may 
well be associated with the system. The radio observations of 
Huang et al. (1994) do show that some emission has diffused 
beyond the optical edge of the galaxy within this region, per- 
haps indicative of previous supernova outbursts. This would 
help to explain the existence of hot X-ray emitting gas also. 
The number of counts within this region make it difficult to 
constrain its spectral properties. 

The 'diffuse' spectrum of Arp 270 as a whole appears 
severely contaminated, both by unresolved point sources, 
and by 'stumps' of source emission left in the source- 
subtraction procedure (see Section 3.1). This is reflected in 
the rather wide confidence errors quoted on the best fit val- 
ues in Table || (remember, no goodness of fit can be quoted 
as a maximum likelihood technique had to be used because 



of the small number of counts). As discussed in Section 3.1 



we attempted to fit a two-component model to the diffuse 
spectrum, comprising a Raymond & Smith hot plasma (rep- 
resenting the hot gas), and a highly absorbed hot (lOkeV) 
component (representing the source contamination compo- 
nent). This two-component fit suggests, as shown in Table || 
that approximately two-thirds of the 'diffuse' emission may 
be truly diffuse gas at a very low (< 0.2 keV) temperature, 
while the remaining third of the emission can be accounted 
for in terms of highly absorbed, hard (10 keV) sources. An F- 
test shows that the improvement in fit quality to the 'diffuse' 
spectrum from including the second component is significant 
at over 99% confidence. 



4.2 Arp 242 

One of the original systems presented by Toomre & Toomre 
(1972) as a classic example of a pair of galaxies undergo- 
ing tidal interaction is Arp 242 (the Mice). Also known as 
NGC 4676, this system lies second in the proposed evolu- 
tionary sequences of both Toomre (1977) and Hibbard & 
van Gorkom (1996), as well as in the present work. Recent 
kinematical work (e.g. Mihos, Bothun & Richstone 1993) has 
strengthened further the tidal interpretation of the Mice's 
structure, and it is now generally agreed that both galax- 
ies, in the throes of a prograde encounter, have their north- 
ern edges moving away from us, such that NGC 4676a's 
tail is on the very furthest side, swinging away from us, 
and NGC 4676b is rotating clockwise, its northeastern por- 
tion, the closest to us (Hibbard & van Gorkom 1996). Hib- 
bard & van Gorkom (1996) find these kinematics in general 
agreement with their Hi velocity data. What discrepancies 
they do see allow them to derive a mass ratio of 2:1 for 
NGC 4676b:NGC 4676a. 

Both the northern galaxy (NGC 4676a) and the south- 
ern galaxy (NGC 4676b) appear to have shapes and colours 
consistent with those of early-type spirals, though the disc 
regions are strongly distorted or absent. Both tails, although 
bluer than the galaxies' central colours, are in agreement 
with the colours of outer-disc regions (Schombert, Wallin & 
Struck-Marcell 1990). They account for 16% of the total Ha 
emission, are quite luminous, containing one-third of the to- 



tal R-band luminosity of the system, and have a high atomic 
gas content (Hibbard & van Gorkom 1996). 

The northern galaxy appears to exhibit a 6.6h~ kpc 
plume of Ha along its minor axis, and the southern galaxy 
possesses an ionized gas bar, as produced in Barnes & Hern- 
quist's (1991, 1996) merger simulations, offset with respect 
to the stellar bar (Hibbard & van Gorkom 1996). Angular 
momentum transfer between the two bars is able to force 
large amounts of gas towards the galactic centre. 

Arp 242 has never been observed in X-rays previously. 
The X-ray image (Fig []) shows rather an amorphous X- 
ray structure. The one source detected (labelled 'A') lies 
further than the ROSAT rms attitude uncertainty (~ 6"; 
Hasinger et al. 1992) from the nuclei of either galaxy, lying 
approximately 40" south of NGC 4676a and 20" west of 
NGC 4676b. It may be, that the source is actually associated 
with the contact region between the two galaxies, and is 
either collisionally heated gas or possibly due to enhanced 
star formation taking place in the dense molecular clouds in 
the contact region. Unfortunately, due to lack of counts and 
limitations on resolution, very little can be said about the 
spectral properties of the source, or its true position. 

Leaving aside this source for the moment, the remainder 
of the emission does seem, in general, to follow the optical 
'heads' of the galaxies, running from the south-east to the 
north-west. The little 'spur' to the south-west, though in- 
teresting, is probably spurious, or if not, unlikely to be asso- 
ciated with the system. The feature to the north-west how- 
ever, may well be associated with the NGC4676a Ha plume 
(Hibbard & van Gorkom 1996) extending in the same direc- 
tion. The correlation of Ha and X-ray features extending 
along the minor axis of galaxies, is often a good indicator of 
a starburst-driven galactic wind taking place, both the Ha 
and the X-ray emission being due to clouds (of different den- 
sities) being shocked by the hot, fast wind from the galactic 
nucleus (Heckman, Lehnert & Armus 1993) . While the fitted 
temperature to the diffuse emission spectrum is in very good 
agreement with fitted temperatures of more nearby, known, 
starburst winds (Heckman 1993; Read et al.1997), the small 
number of counts, make it impossible to be confident of this 
result. 



4.3 NGC 4038/9 

The Antennae, NGC 4038/9 (also Arp 244) is perhaps the 
classic example of a pair of galaxies in gravitational inter- 
action, with spectacular tails spanning nearly 150 kpc. It 
has been the study of many dynamical models (Toomre & 
Toomre 1972; Barnes 1988), and the basic validity of these 
models has been confirmed by studies in neutral (van der 
Hulst 1979) and ionised (Amram et al. 1992) hydrogen gas. 
Ha emission knots (Rubin, Ford & Dodorico 1971), coinci- 
dent with powerful radio features (Hummel & Van der Hulst 
1986; Wolstencroft 1988) are found throughout the central 
parts of both galaxies, and Van der Hulst (1979) showed 
that about 70% of the H I in the system lies in the optical 
tails. Both infrared (Bushhouse & Werner 1990; Vigroux et 
al. 1996) and CO (Stanford et al. 1990) observations indi- 
cate that both the nuclei, and the contact region between 
the two galaxies, are likely to be sites of very active star 
formation. 

The results of the Einstein observation of the Anten- 
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Table 3. The results of the PSPC spectral analysis of Arp 270. Tabulated luminosities are those escaping from the system. Only the 
best model fit is shown. Errors on the spectral fit parameters are lcr for one interesting parameter. 
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Table 4. The results of the PSPC spectral analysis of Arp 242. Tabulated luminosities are those escaping from the system. Only the 
best model fit is shown. Errors on the spectral fit parameters are lcr for one interesting parameter. 



nae (Fabbiano & Trinchieri 1983) were rather inconclusive, 
though the emission was seen to be extended, and contained 
a soft X-ray contribution. Harder emission and a possible 
hard point source at the contact region were also seen. The 
ROSAT PSPC observation sheds far more light on to the 
matter, and full details of the observation, the analysis and 
the results, are given in Read et al. (1995), though a sum- 
mary is presented below. 

The ROSAT PSPC data were processed initially in a 
manner very similar to the method presented here, though 
further, more sophisticated techniques were used in the lat- 
ter stages of the analysis (see Read et al. 1995) . Using a max- 
imum entropy reconstruction technique (Gull 1989; Skilling 
1989), a number of discrete components to the X-ray emis- 
sion were identified (A— G; see Table |B| and Fig. ^). Many of 
these have counterparts at other wavelengths; B and D are 
associated with the northern (NGC 4038) and the south- 
ern (NGC 4039) nuclei respectively. Their properties, the 
facts that they are best fitted with low-temperature and rel- 
atively low-column plasma models, that there are offsets be- 
tween the X-ray positions and the radio positions, and that 
it is known that there is a great deal of supernova activity 
taking place at these sites, suggest that they are bubbles 
of hot extranuclear gas, rather than the central starbursts 
themselves. Features A, C, and possibly E have massive Ha 
knots and radio features as counterparts, and are likely to 
be giant H 11 regions. 



Nearly half of the Antennae's, X-ray flux appears dif- 
fuse, and a two-component spectral fit (see Table |B|), indi- 
cates that almost all of it is genuine, hot (4 x 10 6 K), low- 
metallicity gas (the rest being due to unresolved sources), 
with a total mass perhaps exceeding 10 9 Mq . The bulk of 
this emission appears to envelope the entire optical sys- 
tem, except for the tidal tails. An enhancement (accom- 
panied by a rise in temperature to over 10 7 K) is seen in 
this diffuse emission at the contact point between the two 
galaxies, where it is believed the two discs are colliding. 
This may be due to gas heated by the collision, or to trig- 
gered star-formation in the vicinity. Of great interest are 
the two elongated structures extending to the northwest and 
to the southwest, extending from the galactic discs to radii 
of ~ 30kpc, and culminating in symmetrically-positioned, 
apparently point-like, sources (P and Q). These streamers 
most likely indicate the existence of galactic winds, their 
normal bipolar structure possibly having been disturbed by 



the rapid dynamical evolution of the system (see Section 5.3 
for further discussion of this). Containing « 10% of the total 
hot gas present, these streamers could represent a significant 
loss in mass and energy from the system. 

The point sources at the ends of the two arms are a 
real puzzle. The sources could be related to the arms if they 
were massive (~ 10 8 Mq) objects, ejected from the galactic 
nuclei (see Saslaw, Valtonen & Aarseth 1974), the ejection 
process giving rise to the juxtaposed 'wakes' of hot gas. How- 
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Figure 3. Contours of (0.1— 2.3keV) X-ray emission shown su- 
perimposed on an optical image of Arp 242. Contour levels in- 
crease by factors of two from 7.2 X 10 — 4 ct s" 1 arcmin - 2 .) 



Figure 4. Contours of (0.1— 2.3keV) X-ray emission shown su- 
perimposed on an optical image of Arp 4038. Contour levels in- 
crease by factors of two from 7.2 X 10 — 4 ct s — 1 arcmin~ 2 .) 



ever there are various difficulties with this idea (Read et al. 
1995). Alternatively they could be foreground or background 
objects, for example quasars or stars, completely unrelated 
to the diffuse arms. This latter possibility is supported by re- 
cent optical observations of source P, the optical feature co- 
incident with the source at the end of the northern streamer, 
which indicate that it is a distant (z — 0.155) galaxy with 
an active nucleus (L.R.Jones, private communication). 

The Antennae has also been observed with the ROSAT 
HRI (Fabbiano, Schweizer & Mackie 1996), and the results 
are generally consistent with the PSPC results. A complex 
X-ray structure is seen involving filamentary regions, closely 
following the Ha distribution, emission peaks coincident 
with Hn regions, super-Eddington sources, and prominent 
nuclear sources. There is also some evidence of nuclear out- 
flows and superbubbles. 

4.4 NGC 520 

The nature of the peculiar system NGC 520 was once a puz- 
zle; is it one disturbed galaxy or two interacting galaxies? 
Recent studies have clarified the interacting galaxies inter- 
pretation. Stanford and Balcells (1990) have detected two 
galactic nuclei, just visible in the optical but more clearly in 
their K-band image. The less massive component (by per- 
haps more than an order of magnitude) is the northwestern 
knot, which is optically brighter than the main component. 
The main component is optically weak because it is seen 
edge-on, and the light from its central region is absorbed by 
interstellar dust in its disc (visible as a dark lane in Figure 1 
of Bernlohr 1993). Furthermore, Two hypotheses, either that 
the nearby dwarf galaxy UGC 957 might be primarily re- 
sponsible for the disturbed morphology of a single galaxy 
in NGC 520, or that two interacting disc systems formed 
NGC 520, were tested with numerical simulations (Stanford 
and Balcells 1991). The simulations indicate that NGC 520 



contains two interacting discs which collided ~ 3 x 10 years 
ago and UGC 957 was only involved in the producing of the 
northern half of a tidal tail. 

Tovmasyan and Sramek (1976) found that the compact 
radio source in NGC 520 is situated in the dark lane between 
the two visible parts of the system. Condon et al. (1982) 
later resolved this into a 6" east-west extension, consistent 
with an edge-on disc, lying 3° off east-west. Millimeter-wave 
interferometer maps of the 2.6/xm CO emission (Sanders et 
al. 1988a) show a strong peak at the position of this ra- 
dio source; approximately 1.9 x 1O 9 M0 of molecular gas is 
concentrated in a region approximately 0.8 kpc in size. 

Much of the extranuclear regions of both galaxies within 
NGC 520 experienced a period of enhanced star formation ~ 
3 x 10 8 years ago. The main sequence remnants of this burst 
are the A stars whose features are evident in the optical 
spectra (Stanford 1991). The burst within the less massive 
(the north-western) nucleus occurred slightly later than in 
the extranuclear regions, but the star formation rate has 
returned to a nominal level. The more massive, optically 
hidden nucleus (to the south-east) produces stars at a rate 
of ~ 0.7Mq yr~ and is the current dominant source of star 
formation in this system. The star formation rate within 
this region is ~ 35 times higher than for an isolated disc 
galaxy. This region dominates the mid-infrared flux of the 
system, and probably produces most of the far-infrared flux 
seen in NGC 520. It is still unclear, though, whether the 
massive star component is sufficient to power the activity at 
the very centre, even though multi-wavelength observations 
argue against the presence of a nonstellar compact power 
source within the south-eastern nucleus (Stanford 1991). 

The NGC 520 PSPC data were processed exactly as 
described in section []. Figure ^ shows contours of (0.1—2.3) 
X-ray emission, superimposed on an optical image. Table M 
shows the results of the spectral fitting to the PSPC data as 
described in Section 51. The one source detected in the vicin- 
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Table 5. The results of the PSPC spectral analysis of NGC 4038/9. Tabulated luminosities are those escaping from the system. Only 
the best model fit is shown. Errors on the spectral fit parameters are la for one interesting parameter. 



Figure 5. Contours of 0.1— 2.3 keV X-ray emission from NGC 520 
overlayed on an optical image. The contour levels increase by 
factors of two from 7.2 X 10 — 4 ct s _1 arcmin -2 . 



ity of NGC 520 at a = 01/i24m34.76s, 5 = +03d47m39.7s, 
lies within ~ 5" of the radio source resolved by Condon et 
al. (1982), once the residual error of the ROSAT attitude 
solution (~ 6"; Hasinger et al. 1992) is taken into account. 

It appears that NGC 520 is a very compact X-ray source 
with little in the way of diffuse emission. The Einstein im- 
age (Fabbiano et al. 1992) appears to show the same thing 
though to a much coarser resolution. This compactness can 
be seen more clearly in figure [], a comparison of the ra- 
dial profile of the emission from NGC 520 (crosses) with the 
ROSAT PSPC PSF (line), where very little emission is seen 
beyond a radius of 0.6' (the 95% enclosed energy radius of 



Figure 6. A comparison of the radial X-ray emission from 
NGC 520 (crosses) with the iJOS^T PSPC PSF (line). 

the PSPC PSF). The X-ray emission from NGC 520 is al- 
most consistent with point source emission, and what diffuse 
emission exists, only makes up a very small fraction of the 
total. 

The diffuse emission that could be extracted (using the 
method described in section |^) amounts to only about 20 
counts, and little can be said about its spectral properties, 
except that it appears soft, perhaps suggestive of a very low- 
level, cool corona surrounding the system. This apparent 
lack of diffuse emission in NGC 520 is discussed in detail 
later in Section |5.2j . 

The one point source detected, coincident with both 
the compact radio source resolved by Condon et al. (1982), 
and with the more massive of the two nuclei, as visible in 
the K-band image of Stanford & Balcells (1990), appears to 
dominate the X-ray emission. The power-law fit and the lack 



© 0000 RAS, MNRAS 000, 



12 A.M. Read & T. J. Ponman 



Source Spectrum Analysed 


Net 


log L x 




Best Spectral Fit Results 




R.A. 
(2000.0) 


Dec 
(2000.0) 


Counts 


erg s — 1 
(0.1-2.0 kcV) 


column 
10 20 cm" 2 


photon 
index 


temp. 
keV 


Z 

(solar) 


X 2 (n.d.o.f) 


Integrated 




100.2±10.5 


39.94 


co 7 +47.7 
' -30.5 








22.53 (19) 


Diffuse 




19.4+10.1 


39.12 


00 +83 1 




43+ 21 
u - 4,3 -o.io 


10.01™ 




Source A 01 24 34.76 


+03 47 39.7 


84.7+9.8 


39.92 


9 27+ 11 ' 1 


1 90+ 080 
i ' yu -0.56 









Diffuse + sources 39.99 



Table 6. The results of the PSPC spectral analysis of NGC 520. Tabulated luminosities are those escaping from the system. Only the 
best model fit is shown. Errors on the spectral fit parameters are la for one interesting parameter. 
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Table 7. Sources within NGC 520 detected with the HRI 
minosities are calculated as described in Section 



3.2 



Lu- 
using a 

0.3 keV (and a 3.0 keV) thermal bremsstrahlung model, corrected 
for Galactic absorption. 



of diffuse emission perhaps point to an AGN as the primary 
nuclear pow er source, but NGC 520 is, as shall be seen in 
Section 5.2, far too X-ray dim. No source is detected at the 



position of the north-western nucleus, though an elongation 
of the X-ray contours in this direction is tentatively visible 
in figure |^. 

In contrast to the PSPC image, the HRI image (Fig. 0) 
shows a great deal of detail. Table ^ gives the posi- 
tions, net counts and X-ray luminosities (assuming 0.3 keV 
(and 3.0 keV) bremsstrahlung models) of the detected HRI 
sources. Three sources are detected, the most northerly, HI, 
being centred less than 5" (i.e. less than the residual error in 
the ROSAT attitude solution) east of the secondary (opti- 
cally brighter) nucleus. This source, undoubtedly associated 
with the northern nucleus, appears to be the hardest of the 
sources, as it is the only source detected in the hard band 
(channels 6—11; see Section fo.2| ). The suggested extension 
to the east is interesting, though its existence is tentative. 
Sources H2 and H3 both appear to be much softer than 
source HI, and they follow the bright band of optical emis- 
sion down the north-eastern side of the system. The primary, 
optically fainter, southern nucleus actually lies slightly south 
of H2, and west of H3, where some X-ray enhancement is 
visible. 



4.5 Arp 220 

Arp 220 (also IC 4553, UGC 9913) is the prototypical ul- 
traluminous galaxy, and is the most luminous object in the 
local universe. Being the closest and brightest of these IRAS- 
discovered (Sanders et al. 1988b) systems, it has been exten- 
sively studied across the whole range of the electromagnetic 



Figure 7. Contours of ROSAT HRI X-ray emission from 
NGC 520 overlayed on an optical image. The contour levels in- 
crease by factors of two from 2.2 X 10 — 3 ct s" 



spectrum. The origin of Arp 220's intense power however, 
remains still somewhat of a mystery. On the one hand, the 
huge quantities of dense, warm molecular gas (Scoville & 
Soifer 1991), together with the strong stellar CO absorp- 
tion lines observable in the IR (Armus et al. 1995), point 
to a massive burst of star-formation as the primary energy 
source. On the other, there appears to be a lot of indirect 
evidence for the presence of a quasar obscured by a layer 
of dense gas and dust. Lonsdale et al. (1994), for instance, 
report a OH megamaser in Arp 220 with an OH line peak 
originating in a structure <1 pc across, and they conclude 
that most of the emission from Arp 220 arises from a very 
small region, perhaps on a scale of less than 10 pc, hence 
the quasar nucleus hypothesis. Prestwich, Joseph & Wright 
(1994) suggest that, while the IR spectrum of Arp 220 is 
well fit by a starburst model (Rowan- Robinson & Crawford 
1989), its relative deficit of ionizing photons may indicate 
that it is an evolved starburst, now dominated by an active 
nucleus. They also point out though, that extinction plus 
dust absorption of the Lya photons could result in a much 
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decreased recombination line flux. A third hypothesis for the 
production of such large amounts of power, as discussed be- 
low, involves the idea of a galaxy-galaxy collision. There is 
much evidence for Arp 220 being a merger of two galaxies. 
Majewski et al. 1993, for instance, have observed the very 
close, double nucleus in the near-IR. The optical morphology 
is very similar to that seen in dynamical simulations of late- 
stage mergers of two disc galaxies (e.g. Mihos & Hernquist 
1994). Furthermore, Wright et al. (1990), having obtained K- 
band images of several merging systems, discovered that the 
radial profile of Arp 220 appeared to follow an r 1//4 law, typ- 
ical of elliptical-like, and therefore possibly merger-remnant, 
systems. Within this merger scenario, it has been suggested 
(Harwitt et al. 1987), that the FIR emission and the soft 
X-ray/EUV radiation could be produced by dust heated by 
the collision, and by the cloud-cloud collisions that occur 
when two-gas-rich systems encounter one another. 

Arp 220 has been observed in X-rays prior to ROSAT, 
with the Einstein IPC (Eales & Arnaud 1988). An X-ray 
sources ~ 1' from the position of Arp 220 was discovered, 
and was assumed to be associated with the system. Because 
of the size of the source (thought to be about twice the 
optical size) and the softness of the spectrum, the authors 
concluded that the bulk of the emission could not be coming 
from any obscured active nucleus at the galaxy's centre. 

A thorough analysis of both the ROSATPSPC and HRI 
data from Arp 220 has already been published (Heckman et 
al. 1996), and a discussion of our results in conjunction with 
theirs, is presented here. The ROSAT image (Fig.js]) shows 
quite clearly why the Einstein image was so confusing. Two 
distinct regions of X-ray emission are seen around Arp 220. 
One is centred on the system itself, and appears elongated in 
an approximately east-west direction. The other feature is 
a large (~ 3') amorphous structure, and appears to stretch 
out from Arp 220 to the south and south-west. Two sources 
are detected. Source A lies within the large southern feature 
and source B lies at the centre of Arp 220. 

The northern feature, the approximately east- west 
structure centred on source B, at the centre of Arp 220, is 
large, « 40 kpc long at the distance of 76 Mpc assumed here. 
Interestingly, Heckman et al.'s (1996) HRI image also shows 
an elongation (~ 11 kpc long), though in a more southeast- 
northwest orientation, misaligned by about 30° with the 
PSPC feature. The HRI image also shows, at the very cen- 
tre of the emission, what appears to be a pair of pointlike 
sources of equal brightness, separated by w 5", lying in a 
southeast-northwest orientation, similar to what is seen in 
optical R-band and Ha images. It seems very likely that this 
double structure is related to the prominent dust lane run- 
ning through the centre of the galaxy. The X-ray structure 
is probably an elongated feature which appears as a dou- 
ble source because of the dust lane running perpendicularly 
through its centre. 

The northern east-west elongated structure appears 
very similar to the Ha emission-line nebula visible around 
Arp 220 (Heckman, Armus & Miley 1987; 1990), which, as 
shown by Heckman et al. (1990), appears to be a pair of 
expanding bubble-like structures. As noted by Heckman et 
al. (1996), both the X-ray and the optical features show a 
central maximum, surrounded by a relatively bright, 0.5' in- 
ner region, elongated in the southeast-northwest direction. 
Outside this, a larger structure is visible in both the X-ray 



and the optical, elongated more in the east-west direction. 
Heckman et al. (1996) conclude that both the X-ray and 
the Ha features are the same size (to within 10%). They 
also suggest that both the X-ray and the optical emission 
arise as the result of a bipolar wind, driven out from the 
nucleus over a timescale of a few times 10 7 yr by an ultra- 
luminous starburst (though it is also possible that the wind 
could be driven by a hidden QSO). Heckman et al. (1996) 
also suggest that the ~ 30° misalignment of the outer struc- 
ture may reflect a change in orientation of the system, as 
the encounter has progressed (see Section 5.3.2 ). 

In the previous analyses, the X-ray features seen have 
not been much larger than the optical sizes of the systems 
in question (except perhaps in the case of the Antennae, 
where a more rigorous analysis has been performed; see Sec- 
tion 4.3). In Arp 220 however, the X-ray emission extends 
well beyond the optical confines of the system, and so the 
spectral extraction procedure has been performed in a differ- 
ent way to that described in Section ^. Firstly, a 'northern' 
spectrum was formed, in an identical manner to that de- 
scribed in Section |^, from a circular area slightly greater 
than the optical diameter (actually centred on source B 
and extending out to a radius of 65"). Source spectra for 
sources A and B have been formed exactly as described 
in Section ^[ The 'diffuse' spectrum however, has been ex- 
tracted from a much larger area, chosen to include the south- 
ern emission - a circle of 3' radius, centred halfway between 
sources A and B. The emission from sources A and B has 
been removed, as described in Section ^| 

The 'northern' spectrum and its resultant fits (see Ta- 
ble |8]), therefore refer to the whole of the northern feature, 
the east-west feature surrounding (and including) source B. 
Heckman et al. (1996), in their fitting of the spectrum ob- 
tained from this area, froze the absorbing hydrogen col- 
umn at two values; the column out of our own galaxy 
(3.7x 10 20 cm -2 ), and the maximum column allowable by 
the Hi data of Baan et al. (1987) (5x 10 21 cm" 2 ). Neither 
fit can be excluded. The best fit obtained in our present 
work (though relatively bad) is consistent with the high- 
column fit of Heckman et al. (1996). Consequently, our de- 
rived L x (0.1-2.0keV) of 1.5 x 10 41 erg s -1 is towards the 
higher side of the Heckman et al. (1996) range. Little can be 
said about the different spectral components involved in the 
northern feature except that the central peak (source B), 
is best fitted with a very high-column plasma, and may be 
spectrally harder than the total 'northern' emission, indi- 
cating that the outer, non-source B emission, - the 'bubble' 
emission, is likely to be spectrally, quite soft, as is seen in the 
halos of many starburst and merging systems (e.g. Heckman 
1993; Read et al. 1997). 

The southern feature is very large, about 3' (65 kpc) 
in diameter, bright, accounting for about half of the total 
counts, and appears to be spectrally, quite soft - the best 
fit to the 'diffuse' spectrum (i.e. the entire emission minus 
the emission from sources A and B) (Table |^) being consis- 
tent with a low temperature (6 million K) plasma, absorbed 
merely by the column within our own galaxy. 

Heckman et al. (1996) suggest that this southern fea- 
ture is associated with a background group/poor cluster of 
galaxies, though they cannot exclude the possibility that 
it is in fact, associated with Arp 220. From the apparent 
magnitudes of the galaxies seen at this position in optical 
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survey plates, a redshift for the putative group can be esti- 
mated, and the implied luminosity, size and temperature of 
the southern feature turn out to be quite reasonable when 
compared to known groups or poor clusters. Recent spectro- 
scopic observations of two of these galaxies do show them to 
have similar redshifts [z ~ 0.09; L.R.Jones; private commu- 
nication) comparable to the values estimated by Heckman et 
al. (1996). This shows that there is at least a pair of galax- 
ies behind the southern feature, and strengthens the group 
interpretation. 

In order to investigate the spectral characteristics of 
this southern feature, a 'southern' spectrum was extracted 
from a 1' radius circle centred at the position a = 
15 fe 34 m 53.7 s ,<5 = +23 d 28 m 32 s . X-ray emission coincident 
with the optical galaxies was excluded from the data to 
avoid possible contamination of the spectrum from sources 
associated with the galaxies themselves, and then standard 
spectral models were fitted. The model that best fits this 
southern spectrum is a Raymond & Smith (1977) hot plasma 
model with an absorbing column (5.87±2.74 x 10 20 cm -2 ) 
consistent with that out of our Galaxy, a low temperature 
(0.5 keV), and a very low (< 0.05 solar) metallicity. Fur- 
ther investigation of the temperature-metallicity parameter 
space however, reveals that the temperature is not very well 
constrained; the fit is only marginally worse when the tem- 
perature is ~1 keV, and the metallicity is a few tenths solar. 
At the distance of Arp 220, the best fit model, once scaled 
to allow for the flux lost in the removal of emission associ- 
ated with the background galaxies, gives rise to an emitted 
(0.1-2.0 keV) X-ray luminosity of 1.4 x 10 41 erg s" 1 . The 
equivalent luminosity at a distance corresponding to a red- 
shift of 0.09 (360 Mpc) is 2.9 x 10 42 erg s _1 . 

Assuming first that the southern feature is associated 
with Arp 220, the total X-ray (0.1-2.0 keV) luminosity is 
given by the sum of the luminosities inferred from the spec- 
tral fits to each of the 'southern' and the 'northern' spectra; 
logLx = 41.47 (the same figure is reached by the summa- 
tion of the inferred luminosities of the 'diffuse' spectra and 
sources A and B). If, however, the southern source is not 
associated with Arp 220, then purely the 'northern' spec- 
trum gives \ogLx = 41.18. The true nature of this feature 



is discussed later in Section 5.3 



4.6 NGC 2623 

NGC 2623 (Arp 243), eighth in the Toomre (1977) sequence 
of merging systems, exhibits two large-scale, well defined 
tidal tails, blue colours and distinct regions of star for- 
mation (Schombert et al. 1990). Both Toomre (1977) and 
Joseph and Wright (1985) infer from this morphology that 
NGC 2623 is the product of a tidal encounter between two 
similarly massive disc galaxies. 

Whereas three distinct condensations can be identified 
in optical images of NGC 2623's central region, near-infrared 
images (e.g. Joy and Harvey 1987, Stanford and Bushouse 
1991) reveal a single symmetric nucleus lying between the 
two northern optical features. Joy and Harvey (1987) con- 
clude that the southern condensation cannot be a remnant 
nucleus, and may simply be a counterarm generated by the 
tidal interaction (c.f. Toomre and Toomre 1972). Stanford 
and Bushouse (1991) find that their K-band surface bright- 
ness profile is, as seen in Arp 220, well fitted by an r 1 ' 4 law, 



Figure 8. Contours of 0.1— 2.3keV X-ray emission from Arp220 
overlayed on an optical image. The contour levels increase by 
factors of two from 7.2 X 10 _4 ct s _1 arcmin - 2 . 



characteristic of elliptical galaxies, out to 3kpc. This was 
also seen by Wright et al. (1990), who concluded that the 
stellar light profile of NGC 2623 bears a very close resem- 
blance to that of ellipticals. 

A radio source in the nucleus of NGC 2623 was resolved 
at 6cm (Condon 1980), with a size of 170 pc and a radio 
spectral index such that F u oc z/ -0 ' 8 . The size of the radio 
source, the slope of the radio continuum, and the ratio of 
100 fim to 21 cm flux (Joy and Harvey 1987) are consistent 
with synchrotron emission from supernovae produced in a 
burst of star formation (Condon 1980, Condon and Broder- 
ick 1986). The radio observations also appear to be consis- 
tent with a weak non-thermal emission-line source (c.f. Keel 
1986). 

Strong Balmer absorption lines reveal the presence of 
early-type stars in NGC 2623's nuclear region, making up 
between one-third to one-half of the central optical emission 
(Joy and Harvey 1987). Emission lines such as Hcv and N 
ii, in contrast to the early-type star distribution, are found 
to originate only within the very centre of NGC 2623, and 
the high N n/Ha ratio observed indicates that the nuclear 
emission lines are wholly or partially excited by shocks or a 
weak nonthermal LINER source (Keel 1986). Joy and Har- 
vey (1987) conclude that, since nebular emission lines are 
not associated with the extended distribution of young early- 
type stars, most of the highly ionizing O and early B type 
stars have evolved off the main sequence. A similar conclu- 
sion was reached by Larson and Tinsley (1978), who found 
that NGC 2623 exhibits UBV colours of a galaxy in which 
a burst of star formation occurred w 10 s years ago. 

There are several observations indicating that collisions 
between gas clouds are taking place within NGC 2623. Joy 
and Harvey (1987), for instance, have inferred that the 
merger should produce a region of compressed gas and dust 
~ 2" in extent. This is observed within the central regions 
of NGC 2623. Furthermore, observed strong OH absorp- 
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Table 8. The results of the PSPC spectral analysis of Arp 220. Tabulated luminosities are those escaping from the system. Only the 
best model fit is shown. Errors on the spectral fit parameters are lcr for one interesting parameter. 

a The 'northern' and diffuse spectra have been obtained in a different way from that described in Section H. See text for details. 



tion features imply the existence of a dense concentration 
of molecular clouds in the nuclear region (Baan et al. 1985) . 
Also, CO observations (Sanders et al. 1986, Casoli et al. 
1988) indicate that NGC 2623 has 6 x 10 9 M of molecu- 
lar hydrogen gas, highly concentrated towards the radio and 
infrared nucleus, in a region less than 3kpc in size. 

Once the analysis of the system had begun, exactly as 
described in Section ^, it soon became evident that only a 
small number of counts were present. The analysis was there- 
fore restricted to the energy band (0.1—1.2 keV), in order to 
maximise the signal to noise. Figure ^ shows contours of 
X-ray emission from NGC 2623 not only in the 'full' (here, 
0.1-1.2 keV) energy band, but also in soft (0.1-0.5 keV) 
and hard (0.5— 1.2keV) bands. What is very evident from 
figure |)| is that, as seen in Arp 220, two very distinct compo- 
nents comprise the X-ray emission from NGC 2623; a hard, 
compact nuclear feature, and a cool feature, lying outside of 
the optical confines of the system, above the northern edge 
of the western tail. 

Spectral fitting of the integrated, the diffuse and the 
nucear feature's spectra all proved difficult, due to the 
small number of counts. Only very simple models (thermal 
bremsstrahlung and power law with the absorbing column 
frozen at the Galactic value) could be used to fit these spec- 
tra (see Table ^) , and the maximum likelihood statistic had 
to be used in every case. Thermal bremsstrahlung models 
proved the better fit to each spectrum. Although Source A, 
the nuclear feature, was not actually formally detected at 
the 4cr significance level, we have included it here on the 
basis of it having a signal-to-noise ratio greater than 2. Fur- 
thermore, because of the low number of counts, we have 
not given any spectral fitting results regarding source A. 
Source A's spectrum does appear to be hard however (as 
indicated by Figure []) , and a luminosity is given in Table ^| 
assuming a 1.5 keV bremsstrahlung spectrum absorbed by 
the column out of our Galaxy. 

The X-ray emission from NGC 2623 is split into two 
very distinct components. The nuclear feature (source A) 
is compact, and strong, and appears to be spectrally hard. 
It is coincident with published optical, radio, IR and CO 
emission. It is also seen with the HRI (see Fig. |Io| ), though 
only at a very low flux level (3.6 x 10 -4 cts s" 1 ). 

The second feature, lying to the west of the system and 
coincident with seemingly nothing whatsoever, is of much 



Figure 9. Contours of X-ray emission in three different en- 
ergy bands (top: 0.1— 1.2keV, middle: 0.1— 0.5keV, bottom: 
0.5— 1.2keV) are shown superimposed on an optical image of 
NGC 2623. The contour levels increase by factors of two from 
7.2 X 10~ 4 ct s _1 arcmin -2 . 



more interest. Comparison of the radial profile with the 
PSPC point spread function is inconclusive as to whether the 
source is extended or not. If the feature is point-like, then the 
very low temperature could point to a serendipitous quasar 
or white dwarf as possible candidates. Typical quasar spec- 
tra (a — 2.2±0.2; Branduardi-Raymont et al. 1994; Roche et 



© 0000 RAS, MNRAS 000, 



16 A.M. Read & T. J. Ponman 



Source Spectrum Analysed 


Net 


log L x 




Best Spectral Fit Results 


R.A. 
(2000.0) 


Dec 
(2000.0) 


Counts 


erg s — 1 
(0.1-2.0 keV) 


column 
10 20 cm~ 2 


photon temp, 
index keV 


Z x* (n.d.o.f) 
(solar) 


Integrated 




36.0±10.3 


40.93 


3.28(F) 


49+0-28 
u -* y -0.22 




Diffuse 




19. 2± 9.4 


40.57 


3.28(F) 


20+ 10 




Source A° 08 38 23.54 


+25 45 09 


9.6± 3.8 


40.70 


3.28(F) 


1.50(F) 





Diffuse + sources 40.94 



Table 9. The results of the PSPC spectral analysis of NGC 2623. Tabulated luminosities are those escaping from the system. Only the 
best model fit is shown. Errors on the spectral fit parameters are la for one interesting parameter. 

a No spectral analysis results are given for source A due to the low number of counts. The tabulated luminosity assumes a 1.5 keV 
bremsstrahlung spectrum corrected for Galactic absorption (see text). 



al. 1995) however, observed through the Galactic column to 
NGC 2623, can be shown through simulations, to be far too 
hard to correspond to this soft component. Furthermore, as 
only a few hundred white dwarfs were detected in the entire 
ROSAT all-sky survey, it is statistically very unlikely that 
this feature is due to a foreground white dwarf. It is also very 
unlikely that the soft feature is associated with the western 
tail, as it lies well above it, and this tail is both dimmer and 
less blue than its partner (Schombert et al. 1990). The HRI 
image (Fig. [j^) shows something very interesting however 
as regards this western feature - nothing at all. Nothing 
is seen in the HRI data to the west of NGC 2623 on any 
spatial smoothing scale down to a significance level of 2a. 
The lack of any HRI emission is useful, as regards ascertain- 
ing whether this feature is point-like or not. If the source 
were point-like then we should have expected (based on the 
PSPC count rate) w 0.8 ct ksec^" 1 in the HRI data. This 
is not seen; a 2a upper limit to the HRI count rate at the 
position of the PSPC feature gives 0.2 ct ksec -1 . Hence, ei- 
ther the western feature is very variable, or it is extended, 
on a scale substantially larger than the spatial resolution of 
the HRI. We assume below that the western feature has a 
maximum extent of ~ 2', as suggested by Figure ^. 

The fact that this feature is very soft is reminiscent of 
the X-ray outflows (or winds) seen around many starburst 
galaxies (Heckman 1993; Read et al. 1997). If NGC 2623's 
soft feature is an outflow, then it is a rather different outflow 
from the standard starburst type; apart from it being far 
more luminous and cooler, it appears only on one side of 
the system, and it is a great deal larger (~ 30kpc). 

This large cool, diffuse X-ray structure bears much more 
resemblance in fact, to the feature seen in the previous sec- 
tion; the large amorphous structure to the south of Arp 220. 
a system very similar to NGC 2623 (hence their adjacent 
positions in the present sequence), with comparable IR lu- 
minosities, and similar elliptical-like light profiles (Wright et 
al. 1990). The feature in Arp 220, as discussed in the pre- 
vious section, is believed by Heckman et al. (1996) to be 
associated with a galaxy group or small cluster. In the case 
of NGC 2623, however, nothing is seen in the optical. The 



nature of these features is discussed later in Section 5.3 



Figure 10. Contours of ROSAT HRI X-ray emission from 
NGC 2623 overlayed on an optical image. The contour levels in- 
crease by factors of two from 2.2 X 10 — 3 ct s _1 arcmin - 2 . 



4.7 NGC 7252 

In Toomre's (1977) sequence of merging systems, NGC 7252 
(Arp 226) appears last, its structure, combining a single 
relaxed-looking body with long, faint tidal tails, suggestive 
of a possible end-product of the merger process - an ellip- 
tical galaxy. Schweizer (1978) found many features within 
this galaxy, often compared to the 'atoms for peace' sym- 
bol, which conclusively defined NGC 7252 as a late-stage 
merger. These features included two similar, oppositely- 
moving tidal tails, a single relaxed nucleus, chaotic motions 
within the main body, ripples and loops at the edges of 
the main body, and no nearby perturbing systems. Further 
work showed that the light from NGC 7252 follows an r 1//4 
law (e.g. Stanford & Bushouse 1991), typical of ellipticals, 
it exhibits a post-starburst nuclear population (Schweizer 
1982), and it obeys the Faber- Jackson relationship (Lake 
& Dressier 1986), again characteristic of normal elliptical 
galaxies. Many more recent studies have strengthened the 



© 0000 RAS, MNRAS 000, 



The X-ray Evolution of Merging Galaxies 17 



idea that NGC 7252 is a late-stage merger of two gas-rich 
spiral discs (e.g. Schweizer 1990; Hibbard et al. 1994), and 
it is now considered to be the prototype merging galaxy 
remnant. Perhaps the strongest observational support of the 
merger hypothesis comes from the recent VLA observations 
of NGC 7252 (Hibbard et al. 1994), where, while the tidal 
tails were seen to be very abundant (2 x l{PhT 2 Mo) in 
atomic hydrogen, a fact requiring a gas-rich spiral genesis, 
the main body of the system was seen to be virtually empty 
of H I, and appears indistinguishable from a typical elliptical. 
The molecular gas, in contrast to the atomic gas, is tightly 
bound to an inner rotating disc, of radius ~7", the gap be- 
tween the atomic and molecular gas distributions apparently 
filled by diffuse Ha emission. The H I within the tidal tails of 
NGC 7252 is very interesting: the Mhi/Lb ratio is consis- 
tent with that seen in the outer regions of late-type spirals 
(~ 1; Wevers et al. 1986), and the gas and starlight in the 
eastern tail appear to have the same extent, whereas in the 
northwestern tail, the neutral hydrogen appears to extend 
some 40kpc or so beyond the optical emission. From the 
work of Cayatte et al. (1994), this may indicate that the 
progenitor galaxy that gave rise to the eastern tail may be 
of an earlier type than its partner (Hibbard et al. 1994). 

Computer simulations of collisions between two gas-rich 
discs (Borne & Richstone 1991; Mihos et al. 1993) repro- 
duced the essential features of NGC 7252's appearance at 
first, quite adequately, though a retrograde encounter ge- 
ometry (i.e. galaxy rotations in the opposite sense to their 
mutual orbit), and a very tightly bound orbit, such that the 
galaxies were in contact at the start of the simulation, were 
necessary. These simulations were later however, found to 
be inconsistent with the Hi data of Hibbard et al. (1994), 
and only the more recent simulations of Hibbard & Mihos 
(1995) are able to match the observations (and with a much 
more astrophysically satisfying orbital geometry), explain- 
ing the Hi in terms of the material within NGC 7252's tails, 
having reached a turning point in its expansion, falling back 
to smaller radii. It is thought that approximately half of the 
present tail material 10 9 M Q of Hi, and 2 x 10 9 Z/q of 
starlight) will return to the main body over a timescale of 
less than 1 Gyr. 

NGC 7252 has never been observed in X-rays previously, 
and the ROSATPSPC and HRI observations presented here, 
give us a first glimpse at this important system in X-rays. A 
brief analysis of the PSPC data has already been published 
in Hibbard et al. (1994), and their X-ray image appears 
very similar to our Figure [ll] though we have contoured 
the present image to a greater depth. The X-ray source co- 
incident with NGC 7252 appears very centrally concentrated 
and compact, though a couple of tentative extended features 
are seen, apparently connected to the main body, to the east 
and to the northeast. A comparison of the PSPC PSF with 
a radial profile of the X-ray data centred on the peak of the 
emission, suggests that the source is resolved, though only 
at the 2a significance level. 

The best fit to the integrated spectrum indicates that 
the emission from NGC 7252 may be due mainly to quite 
cool (0.43 keV) gas. As discussed in the last paragraph, the 
emission from NGC 7252 may only just be resolved, and 
consequently, the source-subtraction procedure described in 
Section [] had great difficulties in splitting up this emission 
into source and diffuse emission. In fact, only when a very 



Figure 11. Contours of 0.1— 2.3keV X-ray emission from 
NGC 7252 overlayed on an optical image. The contour levels in- 
crease by factors of two from 7.2 X 10 _4 ct s — 1 arcmin - 2 . 



large radius (4') was used, could the fitting of the diffuse 
spectrum settle at sensible values. As can be seen from Fig- 
ure [ll] however, the use of such a large extraction radius 
merely results in the inclusion of emission from features hav- 
ing nothing to do with NGC 7252, notably fro m fea tures 
along the left-hand edge of Figure [ll] (see Section 5.5). 

Consequently, we have assumed the integrated spectral 
fit results to be the best indicator of the true (0.1— 2.0keV) 
X-ray luminosity; 5.3 x 10 40 erg s -1 . The spectrum is best 
fitted by a low temperature plasma, indicative of hot gas 
enveloping the system (the fact that the fitted absorbing 
column appears larger than the Galactic value probably in- 
dicates some contribution from harder, unresolved sources, 
embedded deeper within the galaxy). 

A number of these sources in fact appear in the HRI im- 
age (Figure |l|), where a great deal of structure is seen. Two 
sources are formally detected close to the central body of 
the system (see Table [n]). Source HI is bright, and appears 
to be coincident with the nucleus, lying less than 7" from 
the position given in the RC3. Features are also seen almost 
symmetrically opposed to the southeast and northwest of 
the nucleus. Source H2 is especially interesting, as firstly, it 
appears to be coincident with the base of the eastern exten- 
sion tentatively seen in the PSPC image (Figure [ll]), and 
secondly, it is the only significant feature detected in the 
hard band. Interestingly, the northeastern feature seen in 
the PSPC image appears also to have a tentative counter- 
part visible in the HRI image (to the top-left of Figure |l2|) . 



4.8 AM 1146-270 

AM1146-270 was selected from the Arp-Madore Catalogue 
of Southern Peculiar Galaxies and Associations (Arp and 
Madore 1987). It has a very complex knotty structure within 
its central 10", as indicated by the B-band image of Smith 
and Hintzen (1991), where seven local maxima, none of 
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Table 10. The results of the PSPC spectral analysis of NGC 7252. Tabulated luminosities are those escaping from the system. Only the 

best model fit is shown. Errors on the spectral fit parameters are la for one interesting parameter. 

a The diffuse spectrum has been obtained in a different way from that described in Section H. See text for details. 



AMI 146-270 received only 601 seconds of PSPC time. 
No significant emission was detected from this system, so 
a 2cr upper limit to the count rate was derived at the po- 
sition given in Smith and Hintzen (1991); 4.47 ct ksec - . 
Assuming a 1 keV bremsstrahlung plasma and Galactic ab- 
sorption (5.18 x 10 20 cm -2 ) then an upper limit to its emit- 
ted (0.1— 2.0keV) X-ray luminosity can be derived; Lx < 
3.9 x 10 39 erg s _1 (the equivalent 5keV bremsstrahlung up- 
per limit is 4.1 x 10 39 erg s _1 ). 



Figure 12. Contours of ROSAT HRI X-ray emission from 
NGC 7252 overlayed on an optical image. The contour levels in- 
crease by factors of two from 2.2 X 10 — 3 ct s -1 arcmin -2 . 
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Table 11. Sources within NGC 7252 detected with the HRI. 
Luminosities are calculated as described in Section [j^, using a 
0.3 keV (and a 3.0 keV) thermal bremsstrahlung model, corrected 
for Galactic absorption. 



which seem to be superimposed stars, are seen. A plume 
or fan (13" in length) extends from the western side of the 
galaxy (see Fig. fil). The colours are very blue, similar to 
those of an irregular galaxy, and many blue knots are seen 
surrounding the nucleus (Smith and Hintzen 1991). This 
galaxy was detected by IRAS and appears to be the site 
of much recent star formation. 



5 DISCUSSION THE X-RAY EVOLUTION 
OF MERGING GALAXIES 

The previous subsections have described the X-ray proper- 
ties of eight systems thought to be at different evolutionary 
stages of a merger of two disc galaxies. These systems range 
from 'approaching' mergers, such as Arp 270, and Arp 242, 
where the discs of the two galaxies are still quite distinct, 
and bridges and tails have only just begun to form, through 
to 'contact' mergers, such as the Antennae and NGC 520, 
where, though the galactic nuclei are still quite distinct, the 
galactic discs have begun to interact strongly, through to 
'ultraluminous' mergers, such as Arp 220 and NGC 2623, 
extremely powerful systems on the point of nuclear coales- 
cence, and finally, to 'remnant' mergers, such as NGC 7252 
and AM 1146-270, relaxed, elliptical-like systems, their vio- 
lent star-formation episodes long-since over. 

Although we are assuming here that these eight systems 
represent an evolutionary sequence, one must be careful not 
to assume that they represent eight different evolutionary 
stages of the same system, for several reasons. Firstly, the 
progenitor galaxies are/were all different in detail, in terms 
of masses, mass ratios, gas richness etc. Secondly, the process 
by which each system has come about, i.e. the encounter ge- 
ometries and parameters, will have all been different. Both 
of these facts are bound to affect the evolution of the sys- 
tem. The fact however, that we have used the appearance, 
length and symmetry of each system's tidal tails as one of 
our main selection criteria, means that we are likely to have 
selected broadly similar systems, i.e. prograde, medium im- 
pact parameter, collisions of similarly-sized disc galaxies. 
Forces need to act for a long time for the formation of strong 
tails, and this requires a predominantly prograde encounter 
(Toomre & Toomre 1972). Similarly, the force needs to be 
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strong, i.e. the impact parameter needs to be small (Barnes 
1992), though not so small as to cause the galaxies to merge 
before the tails have had a chance to grow. Lastly, proper, 
symmetrical tail construction requires that the two masses 
should be roughly equal (Toomre & Toomre 1972). 

Therefore, we should really think of the following evo- 
lutionary discussion in conservative terms, not as snapshots 
of a single system at different points in its evolution, more 
in general terms of how systems appear compared to one an- 
other, at 'approaching', 'contact', 'ultraluminous', and 'rem- 
nant' stages. 

5.1 Previous related work 

Before continuing, it is useful to review here more fully rel- 
evant work on the evolution of interacting systems in other 
wavebands. 

Joseph & Wright (1985) studied the infrared proper- 
ties of a sample of merely intermediate-age mergers (i.e. 
covering a smaller evolutionary period than in the present 
study) including three of the systems within our sample, 
namely NGC 520, NGC 2623, and Arp 220. None of the 
systems covered either showed two distinct discs, as in the 
'nearby' mergers, or appeared as relaxed or elliptical-like 
as NGC 7252. Nevertheless, their results are of great im- 
portance; they observed the infrared luminosity to rise and 
fall along their sequence, suggesting that a massive burst of 
star-formation of luminosity ~ 10 12 Lq occurs within merg- 
ing galaxies as the merger progresses. 

Telesco et al. (1988) studied a large sample of 
comparable-sized interacting galaxies from the Arp-Madore 
Catalogue of Southern Peculiar Galaxies and Associations 
(Arp and Madore 1987). They found that the dust colour 
temperature and the interaction strength are both at their 
largest for systems with the smallest galactic separations. 
This result, along with those of Joseph & Wright (1985) 
above, have, as discussed in Section ^, been used in our se- 
lection and chronological ordering criteria (see Table |l|) . 

More recently, perhaps the most important work con- 
cerning the evolution of merging galaxies has been the work 
of Hibbard & van Gorkom (1996), and it is worth summaris- 
ing some of their more important results here, as many of 
them will have a substantial bearing on the discussion of the 
X-ray results in the next section. 

Hibbard & van Gorkom (1996) presented aperture syn- 
thesis observations of the neutral hydrogen distribution in 
a sample of five systems (four of which are taken from 
the Toomre (1977) sequence, and three of which (Arp 242, 
NGC 520 and NGC 7252) appear here). This they supple- 
mented with wide field, narrow band Ha images and deep 
R-band surface photometry. 

They found that, within the inner regions, the starlight 
in the later systems (including NGC 7252) exhibited 
r 1//4 profiles, a natural consequence of violent relaxation 
(Lynden-Bell 1967; van Albada 1982), and obeyed the Faber- 
Jackson relationship, indicating that as these later systems 
can be dated at about 1 Gyr after the merger (Schweizer 
1996), the primary characteristics of ellipticals are created 
rather rapidly compared to the total merging time-scales. 

The ionized hydrogen is seen, in the early systems (in- 
cluding Arp 242), to form minor axis plumes, similar to the 
wind 'blow-outs' of Heckman et al. (1993), and in the later 



systems, to be confined to the inner kpc, presumably asso- 
ciated with a small molecular gas disc. 

The cold gas is especially interesting, and is seen to be 
widely separated in the later systems, with the molecular gas 
concentrated within the remnant bodies and the atomic gas 
predominantly in the outer regions. Simulations have shown 
that the symmetric forces experienced during a merger can 
force at most half of the outer disc material into a tail, the 
rest of the Hi being forced into the inner regions (Hibbard 
& Mihos 1995). As less than a quarter of the total Hi is 
found within these regions however, Hibbard & van Gorkom 
(1966) concluded that most of the centrally forced gas must 
have been converted during the merger to some other form. 

One possibility for what has happened to this gas, is 
that it has been turned into stars, and indeed, the presence 
of an 'E+A' spectrum (Dressier & Gunn 1983) in NGC 7252 
(Fritze-von Alvensleben & Gerhard 1994), and in the almost 
identical NGC 3921 (Schweizer 1996), indicates that a sig- 
nificant amount of the original atomic gas has formed stars. 
Another possibility for the fate of the atomic gas, is that 
it has been converted into molecular gas, though the fact 
that both NGC 7252 and NGC 3921 have below average 
molecular gas contents for their spiral progenitors (Solomon 
& Sage 1988; Young & Knezec 1989), indicates that there 
has been little net conversion of atomic to molecular gas. A 
third possibility for the lost atomic gas, is that it has been 
heated to X-ray temperatures, either through cloud-cloud 
collisions (Harwitt et al. 1987), or through the energy in- 
put from massive stars and supernovae created in one (or, 
as discussed later, more than one) starburst(s). This neatly 
brings us on to the X-ray evolution. 

5.2 The X-ray evolution 

In Table [12] we have gathered together the emitted 
0.1— 2.0keV X-ray luminosities of t he e igh t sa mple systems, 

Note that in 



calculated as explained in Sections 4.1 



to 4.8 



the case of Arp 220, we have given the X-ray luminosity of 
only the northern feature, the southern outflow/group emis- 
sion has been omitted (inclusion of this emission raises Lx 
to 2.9xl0 41 erg s _1 - see Section [h|). 

Also shown in Table |l|, are the statistical errors (ex- 
pressed as percentages) , based on the number of counts (in- 
cluding background), for each source. There is, in addition to 
this statistical error, a second error which arises from uncer- 
tainty in the appropriate model to use. Since the X-ray lu- 
minosities derived here are band limited to the 0.1— 2.0keV 
band, the effects of considering different acceptable mod- 
els are not very large. Comparing the luminosities inferred 
from power-law and hot plasma models in cases where these 
give fits of comparable quality, we find differences of ~1- 
15% (the larger variations applying to the lowest count-rate 
cases) . For a given spectral model, the effects of varying the 
model parameters within the 90% confidence contour results 
in luminosity variations of the same order. As can be seen 
from Table 12| these uncertainties in Lx are similar in size 
to the simple statistical flux errors. 

In order to see how these X-ray luminosities compare 
with the various multiwavelength properties given in Ta- 
ble [l], what is needed is a plot of the various different lu- 
minosity ratios (Lx/Lb, Lx/Lfir etc.) as a function of 
merger age. As discussed earlier, even establishing a reli- 
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System 



Other names 



log L x 
(erg s" 1 ) 
(% error) 



Arp 270 


NGC 3395/6 


40.49 a 


(10.3 %) 


Arp 242 


NGC 4676 


40.76 b 


(35.3 %) 


NGC 4038/9 


Arp 244 


41.03° 


c ( 2.6 %) 


NGC 520 


Arp 157 


39.94 b 


( 9.5 %) 


Arp 220 


UGC 9913 


41.18 M (11.1 %) 


NGC 2623 


Arp 243 


40. 93 6 


(28.6 %) 


NGC 7252 


Arp 226 


40.67 b 


( 7.6 %) 


AM 1146-270 




<39.61 


e (-) 



Table 12. The merging galaxy sample. Emitted (0.1— 2.0keV) X- 
ray luminosities have been calculated as described in Sections |4.l| 
to E^ , by using the spectral fits to various spectra, as follows 
(quoted errors are purely statistical; see text): 
a Summation of the best 2-component fit to the diffuse spectrum 
plus fits to individual source spectra 
b Best fit to the integrated spectrum. 
c See Read et al. (1995). 

d Note for Arp 220, the 'northern' spectrum has been used, i.e. 
emission from the southern component has been ignored (inclu- 
sion of this component increases Lx to 2.9xl0 41 erg s — 1 ). 
e 2<j upper limit (using a 5keV bremsstrahlung model). 



able chronological ordering is fraught with difficulties, and 
dynamical ages for these systems are considerably more un- 
certain. However, we are interested in broad trends through 
the merger process, rather than precise timescales. 

With this in mind, Fig. [l^ shows the temporal variation 
of the X-ray, the far-infrared, the optical and the radio lumi- 
nosities for the merging galaxy sample. Also plotted is the 
variation of the far- infrared colour temperature, Sm/Sioo- 
All values have been scaled to those of the first system in 
the sequence, Arp 270. Shown for comparison, again scaled 
to the values of Arp 270, are, to the left, the mean values of 
Lfir/Lb, Lx/Lb and Lx/Lfir of both the starburst and 
normal spiral galaxy samples from Read et al. (1997), and 
to the right, the sample of elliptical and SO galaxies from 
Mackie & Fabbiano (1997). When compared to the results 
of Fabbiano et al. (1992), we find that the X-ray luminosities 
of both NGC 7252 and AM 1146-270 (upper limit) are very 
low when compared to typical ellipticals. It should be noted 
however, that Lx/Lb for ellipticals is very La-dependent, 
and so, the 'remnant stage' galaxies' X-ray properties should 
really be compared with those of elliptical galaxies of similar 
optical luminosity. It is seen though, that the X-ray luminos- 
ity of NGC 7252 does appear to be quite low when compared 
to similarly optically bright systems (again from Fabbiano 
et al. 1992), and although only an upper limit to the X- 
ray luminosity from AM 1146-270 could be established, this 
limit is still quite restrictive, and its X-ray luminosity also 
appears to be quite low when compared to systems with a 
similar (very low) optical luminosity (note that very few sys- 
tems appear in this region of Lx — Lb parameter space in 
Fabbiano et al.(1992), so this result is not very secure). 

The interaction timescales in Fig. [l^ are certainly only 
approximate, though in the establishing the evolutionary 
positioning, we have been able to use a variety of point- 
ers. Firstly, for convenience, we have placed time zero at 
the point of nuclear coalescence, chosen here to be at 



the position of the brightest system Arp 220. Now, it is 
thought from its tail kinematics and spectral properties 
(Fritze-v. Alvensleben & Gerhard 1994; Schweizer 1990) that 
NGC 7252 is at ~1 Gyr after nuclear merger - hence its posi- 
tioning. AM 1146-270's colours and morphology suggest that 
it is somewhat older. On the younger side of the evolution, 
comparison of the general morphologies of the sample galax- 
ies with the numerical simulations of Mihos & Hernquist 
(1996) suggest ~650 Myr between an Arp 270- type system 
and nuclear coalescence. Finally, NGC 520, whose simula- 
tion lookalikes seem to occur over halfway in, is thought to 
have experienced extranuclear starburst activity, as is ongo- 
ing in the Antennae, and perhaps the Mice, some 300 Myr 
ago (Stanford 1991). 

The evolution appears to start, as one would perhaps 
expect, at an activity level somewhere between that of nor- 
mal and starburst spiral galaxies. The level of activity, 
whether given by L FIR /L B , L X /L B , Lrad/Lb or Seo/Sioo, 
rises during the evolution, to well in excess of that of typical 
starburst galaxies, and then falls again, back to a very simi- 
lar level to that at the start. Bearing in mind that elliptical 
galaxies may well be the end-products of this evolution, it 
is interesting to note that the final Lx /Lb level is quite 
low when compared to that of typical ellipticals, though it 
is within the accepted activity range. Interestingly, there 
is some evidence that the continued evolution of Lx/Lb of 
even later merger remnants (i.e. moving further to the right) 
is seen to rise gradually to a level in line with that of typical 
ellipticals (Mackie & Fabbiano 1997). The X-ray properties 
of the t hese merger remnants are discussed in more detail in 
Section 



5.4 



What may be of critical importance is a comparison of 
the far-infrared evolution (which we may assume bears a 
very close resemblance to the evolution of the current star- 
formation rate, SFR) with the theoretical work of Mihos & 
Hernquist (1996). These authors used numerical simulation- 
sto investigate the gas dynamics and starbursts in mergers 
similar to those considered here (i.e. of systems of compara- 
ble mass) . They found that the actual structure of the galax- 
ies plays a dominant role in the evolution of the SFR. Two 
merging scenarios, identical except for the make-up of the 
systems, were considered. One consisted of a merger of pure 
disc/halo systems, systems completely lacking a bulge com- 
ponent, the other consisted of a merger of disc/bulge/halo 
systems, each dense bulge component having a mass some 
three times that of the disc component. It was found that the 
gaseous inflows (which are related to the SFR, and hence to 
the FIR luminosity) are strongest when galaxies with dense 
central bulges are in the final stages of merging. 

Firstly, let us discuss the SFR evolution of the non- 
bulge mergers. The SFR in these systems is seen to peak 
at 40—60 times field spiral values at a time some 500 Myr 
before nuclear coalescence, this peak lasting for about 150 
Myr, i.e. from a time around the position of Arp 242, till 
soon after the position of the Antennae. As these starbursts 
deplete the nuclear regions of gas, the starbursts die out, and 
the SFR is seen to drop again close to its quiescent value at 
about the position of NGC 520. During nuclear coalescence 
(at time zero), only a relatively weak starburst is produced 
(an increase in the SFR by about an order of magnitude, and 
lasting for around 50-100 Myr). The SFR history of the bulge 
mergers however, is very different and bears a very close 
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Figure 13. The evolution of the merging galaxy sample. Plotted is the temporal variation of the X-ray, the far-infrared, the optical and 
the radio luminosities for the merging galaxy sample. Also shown is the variation of the far-infrared colour temperature, 560/5100. All 
values have been scaled to those of the first system in the sequence, Arp 270. Also shown for comparison, are, to the left, the mean values 
of Lpxn/Lg, Lx/Lg and Lx / Lfir of both the starburst (stars) and normal (diamonds) spiral galaxy samples of Read et al. (1997), 
and to the right, the Mackie & Fabbiano (1997) elliptical and SO sample. 



resemblance to Fig. |l3|. The starbursts that develop after 
the initial encounter are much weaker than those seen in the 
bulgeless case. The SFR is increased by only a factor of a few 
(less than an order of magnitude), though, as the gas is not 
significantly depleted, it lasts longer, from about Arp 270, or 
just after, to just before Arp 220. Consequently, at close to 
time zero, far more gas exists, and the rapid collapse of this 
gas leads to a very violent starburst, with an SFR some 150 
times that of a field spiral, lasting for only 50 Myr. In both 
the bulge and the non-bulge cases, continued infall of the 
tidal debris is unable to fuel any significant star-formation, 
and the remnant after nuclear merger evolves passively at a 
SFR level similar to that of field spirals. 

In reality, galaxies may have a wide range of bulge-to- 
disc ratios, and mergers occur between galaxies of all types. 
Although our merging sample has not been chosen with a 
constant bulge-disc ratio in mind, all the systems seem to 
show the signature of high-bulge progenitors, when com- 
pared to Mihos & Hernquist's (1996) work. The early sys- 



tems (Arp 270, Arp 242 and the Antennae) show only mod- 
est increases in Lfir, perhaps by a factor of ten, above 
normal field spirals. Similarly, the 'ultraluminous' systems' 
progenitor galaxies must have had large bulge-to-disc ratios, 
because of their massive inferred SFRs. This does throw up 
further caveats regarding the difficulty of an evolutionary 
analysis. Firstly a disturbed morphology does not necessar- 
ily mean that a very large burst of star-formation is taking 
place, and secondly, the main requirements for the existence 
of an ultraluminous starburst may be specific encounter ge- 
ometries and/or specific galaxy structures (Hibbard & van 
Gorkom 1996). 

Comparing the X-ray and far-infrared evolution, one 
sees that although both rise and fall along the merger se- 
quence, the X-ray does so to nowhere near the extent of the 
far-infrared (and, to a lesser extent, the radio). This leads 
to the drop by well over an order of magnitude in Lx /Lfir 
in Fig [l3| Before concentrating on this important feature of 
the evolution however, the striking 'glitch' in the Lx/Lb 
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evolution associated with the system NGC 520, should first 
be discussed. 

NGC 520, the 'second brightest very disturbed galaxy 
in the s ky' (Arp 1987), while once thought, as discussed in 
Section 4.4 , to be a single exploding galaxy (e.g. Arp 1967) , 
is now believed, after receiving so much attention over the 
years, to be a merger-induced system (e.g. Stanford & Bal- 
cells 1990, 1991; Bernlohr 1993). Its far-infrared properties, 
as seen in Fig. Il3l appear to sit nicely on the evolutionary 
sequence. Its X-ray properties however, do not. NGC 520 
is very underluminous in the X-ray, with an Lx/Lb ratio 
similar to that of normal (i.e. non-starburst) field spirals. 
Note that NGC 520's very low Lx/Lfir ratio (similar to 
that of Arp 220 and NGC 2623) is very unlikely to be due 
to the same phenomenon as in the cases of these ultralumi- 
nous mergers, as they have high Lx / Lb ratios, NGC 520 
does not. It is, quite simply, X-ray dim. 

This can be seen very clearly in Fig. [li], where the X- 
ray and far-infrared luminosities of the eight merging galax- 
ies are plotted (AM 1146-270's upper limit as an arrow), 
together, for comparison, with Read et al.'s (1997) samples 
of normal and starburst spiral galaxies. Best fit regression 
lines (excluding the AM 1146-270 upper limit) to the merg- 
ing galaxy sample data are plotted both inclusive and ex- 
clusive of NGC 520. Also plotted is the best fit regression 
line to the Read et al. (1997) sample. 

NGC 520 stands out quite obviously in Fig. [13. Al- 
though very far-infrared bright, as bright as the bright- 
est non-interacting starbursts, it is X-ray underluminous by 
nearly an order of magnitude. Why? There is quite a lot of 
evidence to suggest that NGC 520 is a quite different system 
from the rest of this merging sample. 

A strange feature of NGC 520 is the fact that the disc is 
still remarkably intact; a rotating disc of neutral hydrogen 
of mass a few xlO 9 M0 and radius ~ 20kpc is seen within 
the centre of NGC 520 (Hibbard & van Gorkom 1996). Opti- 
cal isophotes indicate that both nuclei are embedded within 
this disc. The disc has survived the merger so far very well. 
What is to say it that it will not survive the nuclear coales- 
cence? Why NGC 520's disc has survived so well is somewhat 
of a mystery, though numerical simulations (e.g. Quinn et 
al.1993) suggest that the mass ratio may be very important. 
For dynamically cold discs to be significantly disturbed, vi- 
olently fluctuating fields such as those accompanying mod- 
erate mass-ratio interactions (Mi jMi > 0.3) are required. 
For NGC 520, whilst molecular linewidth studies (Sanders 
et al. 1988a) and stellar velocity dispersion measurements 
(Stanford & Balcells 1990) suggest mass ratios close to unity, 
rotation curves (Bernlohr 1993) give a mass ratio of 10:1. 

While the presence of a large gaseous disc within the 
merging system requires at least one of the progenitors to 
have been gas rich, Hibbard & van Gorkom (1996) ob- 
serve the extensive plumes within NGC 520 to be gas- 
poor (Mjjj/Lh<,0.1), suggesting that NGC 520 came about 
through the merger of gas-poor progenitor discs (SO— Sa 
galaxies). The most likely scenario offered by Hibbard & 
van Gorkom is that the NGC 520 system is a result of an 
encounter between a gas-rich system with an extensive disc 
and a gas-poor system, such as an SO or Sa galaxy. They also 
conclude though that NGC 520 shows strong evidence that 
mergers may not necessarily destroy the gaseous discs of the 
progenitors and therefore need not evolve into ellipticals. 



Another related possibility, given that extensive-disc (i.e. 
near-bulgeless) systems may be involved, is that NGC 520 
may be the result of a merger between two bulgeless sys- 
tems, i.e. as in the evolution of Mihos & Hernquist (1996) 
discussed earlier, though its relatively low SFR would mean 
that it must be at an evolutionary stage sometime after 
the early strong starburst seen in these simulations. Either 
way, because of the different encounter progenitors involved, 
whether in terms of mass ratio, gas richness or bulge-to-disc 
ratio, NGC 520 does appear to be on a significantly different 
evolutionary path to the rest of the sample. 

Turning now to the X-ray evolution of the sample as a 
whole, it is worth studying Fig. [l^ and Fig. [li] together. A 
rather tight Lx'-Lfir. relationship is seen in Fig. [w| run- 
ning quite smoothly and continuously from quiescent nor- 
mal spirals (diamonds), through starbursts (stars) to very 
active merging galaxies (filled squares). There does appear 
however, to be a gradual flattening of the slope at higher 
activities, and this is borne out by the large change in the 
slopes of the regression fits, from 0.80, for normal spirals, 
to 0.38 (0.31 when NGC 520 is discounted) for the merger 
sample. 

Why does this flattening of the Lx'-Lfir relationship 
occur? Also, why does the large drop in the X-ray to far- 
infrared lumnosity ratio, seen in Fig [l^, occur?. We believe 
the answer to both questions to be the same, as both the 
Lx'-Lfir flattening, and the Lx/Lfir drop are really just 
the same effect manifesting itself in two different ways. As 
the merger progresses, starburst activity is seen to rise and 
fall. The increases in both the far-infrared and the X-ray 
luminosity are both due to the increased starburst activity. 
While the infrared flux is primarily due to dust heated by the 
massive stars within the starburst, the rise in X-ray flux is 
limited, as much of the input energy from the large number 
of supernovae and stellar winds created within the starburst, 
is lost. It is most likely that this energy is lost in the form 
of kinetic energy - kinetic energy associated with the huge 
gaseous ejections seen in the X-ray in the most luminous 
members of the merging galaxy sample. We now move on to 
discuss these diffuse structures. 



5.3 Diffuse structures 

Diffuse structures, as described in the notes on the individ- 
ual systems, are observed throughout almost the whole of 
the evolutionary sequence. Arp 270's two point sources ap- 
pear to be embedded in a small amount of diffuse emission, 
though this emission does not seem to extend very far (if 
at all) beyond the optical confines of the system. The two 
sources are both best-fitted by warm plasma models, how- 
ever, suggestive of hot gas. Arp 242's distance makes it hard 
to make any firm judgements about the extent of any dif- 
fuse structures. The fact that the source appears extended 
however, along with the fact that the emission is signifi- 
cantly softer than in Arp 270, is very suggestive of signif- 
icant amounts of hot gas. Furthermore vigorous disc-wide 
star-formation is seen in Ha within this system, along with 
plumes of gas extending along the minor axis, possibly co- 
incident with the north-west X-ray feature. The coincidence 
of Ha and X-ray features extending along the minor axis of 
this system points towards the existence of a galactic wind 
or 'blowout' (e.g. Heckman 1993), NGC 4038/9, as discussed 
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Figure 14. Plot showing the ROSAT (0.1— 2.0keV) X-ray luminosity versus far-infrared luminosity, Lpm, calculated as described in 
Section |i| for all members of the merging galaxy sample (filled squares) (AM 1146-270's upper limit is shown by the arrow). Also shown, 
for comparison, are the results for the samples of normal (diamonds) and starburst (stars) spiral galaxies from Read et al. (1997). Best 
fit regression lines (excluding the AM 1146-270 upper limit) to the merging galaxy sample data are plotted both inclusive (dashed line) 
and exclusive (solid line) of NGC 520. Also plotted is the best fit regression line to the Read et al. (1997) sample. 



briefly here, and in detail in Read et al. (1995), shows a great 
deal (~ 1O 9 M0) of low temperature (0.4 keV) gas, not only 
enveloping the two discs, in what appears to be a hot corona, 
but also apparently ejected from the system to the north 
and south as a pair of galactic winds. NGC 520 is very com- 
pact, and shows little evidence for any significant amount of 
diffuse gas. However, as we have already discussed, this sys- 
tem may well be an oddity. Arp 220, the most active of the 
sample, shows a 40 kpc long, east- west structure, coincident 
with a 'double-bubble' Ha emission-line nebula (Heckman, 
Armus & Miley 1987;1990). This X-ray/optical emission is 
generally believed to be due to a bipolar wind driven out 
from the nucleus by an ultrapowerful starburst. The large 
soft source to the south of the galaxy is apparently connected 
to Arp 220, though it may in fact be associated instead with 
a background galaxy group, as discussed in Section ^fj. How- 
ever, NGC 2623, a very similar system to Arp 220, exhibits 
a near-identical feature, which seems a remarkable coinci- 
dence. NGC 2623's very soft (~ 0.2 keV) feature appears 
to be truly diffuse, given that it is not seen with the HRI. 
NGC 7252, a near elliptical-like remnant, shows some evi- 
dence for an X-ray halo. This halo of hot gas surrounding 
NGC 7252, however, appears small. And finally, only an up- 
per limit could be calculated for AM 1146-270, due to the 
very short observing time. It is however, quite a restrictive 
limit, and places the activity level of AM 1146-270 close to 
that of NGC 7252, and definitely far less than that of the 
ultraluminous systems. 

The time evolution of the temperature, the size, and the 
Lx(diff)/Ls values of these diffuse structures are summa- 
rized in Fig. [l5| All values have been scaled to those of the 
Antennae since Arp 270's diffuse emission is only suggestive. 
Two sets of points are plotted for Arp 220, one (the higher 
set) for the case where the southern source is assumed to be 



associated with Arp 220, and the other (the lower set), for 
the case where it is not. In a discussion of the trends seen 
within this figure, one must be very conservative. The ex- 
traction of truly diffuse emission and the fitting of models to 
these spectra has been, as described earlier, difficult. Even 
so, there is evidence that, as a merger progresses, hot gas is 
ejected from the system. This gas is seen to grow both in 
size and in X-ray brightness. Curiously, by the late stages 
of the evolution, there is little evidence for much gas. Little 
can be said about the temperature evolution, except that it 
appears to remain low (< 0.5 keV) throughout the merger. 

A major problem with forming ellipticals from spirals is 
to account for the disappearance of the gas from the original 
spirals. Star formation is too inefficient a process to suppose 
that it is all converted into stars, and so the bulk is probably 
blown away by supernova explosions and stellar winds, and 
ultimately, galactic winds. We believe that we are seeing this 
ejected gas in X-rays, notably in the Antennae, Arp 220 and 
NGC 2623. 

A great range of gaseous structures are observed within 
the sequence, and three fundamental (and chronological) 
questions need to be addressed; When and where does this 
gas ejection process start? What is the true nature of the 
huge one-sided structures seen in the ultraluminous sytems? 
Why is there only marginal evidence for a hot halo in the 
remnants? 



5.3.1 Early stages 

Even in a pre-merger system such as the Antennae, a great 
deal of hot gas is seen, some of it being forcibly ejected from 
the system. What about still earlier systems? One would 
perhaps intuitively think that no ejection could occur before 
the galaxies first encounter each other (i.e. before Arp 270), 
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Figure 15. The evolution of the diffuse structures seen within the merging galaxy sample. Plotted are the temporal variations of the 
temperature (stars), the size (in terms of the maximum extent) (squares) and the values of Lx(diff)/Ls (diamonds). All values have 
been scaled to those of the Antennae. Upper limits are given by arrows. Two sets of points are given for Arp 220, one (the higher set) 
for the case where the southern source is assumed to be associated with Arp 220, one (the lower set), for the case where it is not. In the 
case of NGC 520, it proved impossible to extract any believable information regarding any diffuse component, and no values are plotted. 



and this idea is strengthened by the modelling work of Mihos 
& Hernquist (1996) in which no significant star-formation 
occurs before merging galaxies first meet. Arp 270 appears 
to be at such an early stage in its encounter, it just may 
be on the verge of gaseous ejection. Remember that this 
system is at such an early starburst episode that, although 
star-formation is taking place at a level between that of typ- 
ical normal and starburst spirals, it may be occurring well 
away from the nuclei, the gas not having had time yet to 
be compressed into the nuclear regions. Arp 242, at a later 
stage, is almost certainly ejecting hot gas into its halo. Its 
distance makes it difficult for us to say for sure, but its X- 
ray structure, and temperature, along with the coincidences 
of Ha features, point towards 'blowout '-type ejection taking 
place within this system. 



5.3.2 Mid-stages 

These ejections of hot gas are likely to grow as the merger 
progresses, as an increasing amount of gas is compressed into 
the nuclei, giving rise to more and more supernova explo- 
sions. How then do these gaseous structures appear at the 
peak of merger-induced activity? Arp 220 and NGC 2623 
are very similar systems. Both are extremely bright, espe- 
cially in the infrared and in the radio, and are members 
of the class of ultraluminous far-infrared galaxies (FIRGs) 
discovered by IRAS (Sanders et al. 1988b). Both appear 
to be at an evolutionary stage corresponding to the short 
(50— 100 Myr) violent starburst period, seen in the Mihos & 
Hernquist (1996) simulations. At first glance, both appear 
very similar in the X-ray, with a large, soft, gaseous feature 
apparently connected to the galaxy, but situated some way 
from it, and a harder, more compact source centred on the 
galaxy itself. Furthermore, the same phenomenon is seen in 
ROSAT observations of a third violently merging, ultralu- 
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minous FIR galaxy, NGC 6240 (Fricke & Papaderos 1996). 
Again, in this galaxy, which is very similar to both Arp 220 
and NGC 2623 in terms of its extreme luminosity and po- 
sition in the merging evolution, a massive (~ 50kpc), one- 
sided, diffuse feature is seen in the soft band image, whereas 
in the hard band image, only a point-like nuclear source is 
visible. 

The similarities between these three systems strongly 
suggest that the large soft features are associated with the 
mer ging galaxies. On the other hand, as discussed in Sec- 



tion 4.5, there is evidence for a group of galaxies at z ~ 0.09 
lying behind the soft feature in Arp 220. The properties of 
the diffuse X-ray emission from galaxy groups have been es- 
tablished over the past few years as a result of studies with 
ROSAT (Mulchaey et al. 1996, Ponman et al. 1996). The 
temperature of the hot gas is typically ~ 1 keV, and the 
luminosity (for Ho =75 km s - Mpc -1 ) ranges widely from 
10 40 ergs -1 to 10 43 ergs -1 , with a strong positive correlation 



with temperature. Unfortunately, as reported in Section 4.5 



the spectral properties of the diffuse source in Arp 220 are 
not very well constrained. At the best fit temperature of 
~0.5keV, the derived luminosity of 2.9 x 10 42 ergs -1 would 
be very high for a group, but if T ~ 1 keV (as is allowed 
by the data) it would be quite reasonable. Hence the group 
hypothesis for the soft source cannot be ruled out by the 
available data. 

We have a choice, then, between two coincidences. Ei- 
ther the source to the south of Arp 220 is a background 
group, in which case it is a remarkable coincidence that 
very similar features (with no sign of any background galax- 
ies) are seen in the two similar systems NGC 2623 and 
NGC 6240, or the source is associated with Arp 220 it- 
self, in which case (unless we are prepared to invoke non- 
cosmological redshifts) it must be a coincidence that there 
is a group, or at least a pair, of galaxies behind it. 

One might hope, if these features are indeed galactic 
outflows, that one would see some corresponding emission 
from them at other wavelengths. In the case of NGC 2623, 
the soft X-ray source extends well outside the optical galaxy, 
but the Ha+Nn emission-line image (Armus, Heckman and 
Miley 1990) shows very little structure in the direction of 
this outflow. It is believed that both soft X-ray emission 
and optical line emission are produced by the shocking of 
cool clouds embedded in a hot wind, and most normal sys- 
tems where hot outflows are believed to exist (e.g. NGC 253, 
M82) show X-ray and optical line emission features of rea- 
sonably similar size (see Heckman et al. 1993). Arp 220 is 
a case in poin t, its optical emission line halo, as discussed 
in Section |4.5| , is approximately coincident with the smaller, 
east-west X-ray halo. Nothing is seen however, coincident 
with the huge southern extension. 

What are these features? If they truly are hot gas 
ejected from the merger, then several questions arise. Why 
do they have no counterparts at other wavelengths? Why 
are they one-sided? 

The Mihos & Hernquist (1996) simulations suggest that 
in a merger between two similar-sized, gas-rich, large-bulge 
systems (as we believe we have sampled here), two, quite 
different bursts of star-formation occur. The first occurs as 
the sytems first meet each other (i.e. from Arp 270 to some- 
where close to the Arp 220 epoch) , and is of a comparatively 
low level (a level, in fact, similar to that of nearby star- 



bursts, such as NGC 253 and M82), and lasts for a compar- 
atively long time. The second starburst is far more intense, 
and lasts for an extremely short time (~ 50Myr), occurring 
around the time of nuclear coalescence (our time zero), i.e. 
at the Arp 220/NGC 2623 epoch. This starburst is very dif- 
ferent from those experienced by the earlier 'approaching' 
and 'contact' mergers. 

Consider now the possible evolution. The first starburst, 
because of its similarity to M82-types, is likely to produce 
similar bipolar outflows, as are observed, both in the Anten- 
nae (Read et al. 1995), and in Arp 220 (Section 4.5; Heck- 
man et al. 1996). Interestingly, however, in both cases, there 
is evidence that the 'bipolar' nature of these outflows has 
become distorted by the evo lution of the merger (as dis- 

and in Read et al. (1995) and 



cussed in Sections 4.3 



and 4.5 



Heckman et al. 1996)). An alternative explanation for their 
appearance may relate to the creation of 'unipolar' outflows, 
as discussed below. 

Perhaps the strangest facet of the large, soft features 
seen in the most lumunious systems is their unimodality. 
Why do the structures seen in Arp 220, NGC 2623 (and 
NGC 6240) appear only on one side of the system? The 
answer may lie in the fact that these systems are evolv- 
ing rapidly. Mac Low & McCray (1988), using the Kompa- 
neets (1960) (thin-shell) approximation, numerically mod- 
elled the growth of superbubbles: large thin shells of cold 
gas surrounding a hot pressurized interior - the progenitors 
of galactic winds - in various stratified atmospheres. Among 
their findings, they discovered that superbubbles blow out of 
the H I layer when they grow to a radius of between one and 
two scaleheights. The swept up shell then accelerates out- 
wards and fragments, due to Ray leigh- Taylor instabilities, 
and the wind is then able to escape into the inter-galactic 
medium at velocities of several thousand km s -1 (see Heck- 
man et al. 1993). 

However, Mac Low & McCray (1988) also find that 
these bubbles will blow out on one side only of a disc galaxy 
if their centres are more than 50—60 pc from the centre of the 
disc. Now, in a (relatively) dynamically stable system, such 
as M82 or NGC 253, the starburst is found to be very sym- 
metrically positioned with respect to the galactic disc, and 
bipolar structures are seen in the X-ray (Watson, Stanger 
& Griffiths 1984; Fabbiano 1988; Pietsch 1992; Strickland, 
Ponman & Stevens 1997; Read et al. 1997). In a rapidly- 
evolving ultraluminous merging system such as Arp220 or 
NGC 2623 however, the central burst is highly unlikely to be 
so centrally positioned with respect to the quickly-moving 
and highly distorted gaseous components. There will then be 
a single direction of steepest pressure gradient, along which 
the bubble will expand most rapidly, leading to a one-sided 
blowout . 

It is also possible that this might apply to earlier sys- 
tems within the present sample, as alluded to above. The 
distorted 'bipolar' structure seen in the Antennae, for in- 
stance (Section i.c; Read et al. 1995), may not be a bipolar 
wind after all, but rather a pair of unipolar winds, one pro- 
duced from the northern starburst, one from the southern. 
There is some evidence to support this, both in terms of the 
fact that the individual discs have become very distorted, 
and that the individual nuclear starbursts appear to be off- 
set (see Read et al. 1995). The same could be also true in 
the case of the 'double bubble' structure seen in Arp 220. 
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Two individual unipolar winds, as opposed to one bipolar 
wind may have produced this structure. 

Mean physical properties for these hot gas structures 
can be inferred from our spectral fits if we make some as- 
sumptions about the geometries of the emission. We have 
assumed that both the Arp 220 southern feature and the 
NGC 2623 feature are spherical in nature, on the basis of the 
fact that both the Arp 220 and the NGC 6240 feature (Fricke 
& Papaderos 1996) appear quite round (the NGC 2623 fea- 
ture appears quite round also, though we are un sure whether 
it has been resolved with the PSPC (see Section pi~6| )). This is 
the simplest assumption to make and indicates the largest 
volume over which the gas is likely to be distributed. Us- 
ing the volumes derived for this 'bubble' model, the fitted 
emission measure nn\V (where n is the 'filling factor' - the 
fraction of the total volume V which is occupied by the emit- 
ting gas) can be used to infer the mean electron density, n e , 
and hence the total mass Mgas, thermal energy E^ and 
cooling time i coo i of the gas, and also the mass cooling rate 
M coo j and adiabatic expansion timescale iexp- These values 
for the Arp 220 southern feature and the NGC 2623 eastern 
feature are given in Table [l3| along with, for comparison, 
the corresponding values for the Antennae corona (Read et 
al. 1995), and for the NGC 7252, assuming that entire in- 
tegrated emission from NGC 7252 is in the form of a hot 
gaseous halo. 

If, as we are assuming, we are seeing some sort of out- 
flow within these ultraluminous systems, then if n w 1, they 
are presumably bubbles of gas at 2—6 x 10 6 K flowing out 
from the system. If this is the case however, this gas, un- 
less some additional force is confining it, should expand at 
its sound speed (a few hundred km s _1 ), which is too slow 
when compared to the size of these features (w 50kpc)and 
the timescale of the starburst (w 50Myr). Alternatively, 
the soft X-ray emission we see could arise from material, 
e.g. clouds, shock heated by a fast wind (see Heckman et 
al. 1993). Although this hypothesis is supported by recent 
theoretical modelling of starburst winds (e.g. Suckov et al. 
1994), the actual value of the filling factor of the gas emit- 
ting in the ROSAT band, remains unclear, some simulations 
(e.g. Suchkov et al. 1994) predicting a low value, others (e.g. 
D.K.Strickland, private communication) predicting a high 
value. In the low filling factor case, the ROSAT band X-ray 
luminosity of the clouds should far exceed that of the wind 
fluid itself. In the high filling factor case, the wind fluid 
contibutes a very significant fraction to the ROSAT band 
emission. 

In contrast to the corona surrounding the Antennae 
galaxies, therefore, which appears to be bound to the sys- 
tem (n « 1; see Read et al. 1995), the offset features seen 
in Arp 220 and NGC 2623 are (like the Antennae outflows) 
likely to be due to a combination of dense clouds, shocked 
by a hot, fast unbound wind, and the wind fluid itself. Pro- 
vided that this material does escape from the system, then 
the associated mass loss rates can be calculated from the 
inferred gas masses, if we assume a velocity for the mate- 
rial. This gives 140 (for Arp 220) and 50 (for NGC 2623) 
X\/V v wooMq y r_1 , where uiooo is the velocity, in units 
of 1000 km s . Similarly, the rate at which kinetic en- 
ergy is lost can be estimated to be 4.8 (Arp 220) and 1.7 
(NGC 2623) xl0 43 y»7W 1 oooerg s _1 . 

If n were « 1, i.e. if we were seeing the hot wind itself, 



rather than shocked clouds with !)< 1, then this mass- loss 
could remove a substantial amount of gas (~ 10 10 Mq) from 
the merger (a requirement, if these objects are to eventually 
resemble gas-poor ellipticals), even over the short starburst 
timescale of 50 — 100 My. Furthermore, the rate at which 
kinetic energy could be lost turns out to be quite compa- 
rable to the total inferred supernova luminosities of these 
systems (calculated using the method of Read & Ponman 
(1995)); 28 and 6 xl0 43 erg s" 1 , respectively for Arp 220 
and NGC 2623. Hence, if 7] were ~ 1, then almost all of the 
mass and around a quarter of the energy injected into these 
systems by supernovae could escape in the wind. 

If n is substantially less than unity, the above results 
still stand, though they refer to the mass and energy lost 
in the clouds (provided, of course, they escape the system). 
This reduces the mass and energy loss rates calculated, but 
there is of course additional mass and energy loss associated 
with the wind fluid itself, which is invisible to ROSAT. Sim- 
ulations (Suchkov et al. 1994; D.K.Strickland, private com- 
munication) appear to show though, that the lower density 
phases of the gas do not carry a significant fraction of the 
mass. 

It is also interesting to compare the values calculated 
above with the total energy and mass deposition rates, cal- 
culated (as in Heckman et al. 1993) in terms of the starburst 
bolometric luminosity Lboi (which in the cases of Arp 220 
and NGC 2623 we may assume is approximately equal to 
Lfir)- The energy deposition rates E = 8 x 10 42 Lboi,n erg 
s _1 (where Lboi.n is in units of lO 11 !/©), are 6.9 (Arp 220) 
and 1.5 (NGC 2623) xlO 43 erg s _1 , and the mass deposition 
rates M = 3L bo i,ii M Q yr" 1 , are 25.9 (Arp 220) and 5.8 
(NGC 2623) M© yr" 1 . 

As a final point, note that in the above discussion, we 
have only considered the Arp 220 southern feature, not the 
northern 'double-bubble' emission. Inclusion of this emis- 
sion into the analysis, i.e. performing the same calculations 
to t he re sults of fitting the total 'diffuse' spectrum (see Sec- 



tion 4.5), gives rise to estimated gas masses and thermal 
energies almost a factor of two greater than those given in 
Table [l3| Hence if this material is also escaping, it will ap- 
proximately double the mass and energy loss rates discussed 
above. 



5.4 Final stages 

The X-ray properties of the end products of our evolu- 
tion, NGC 7252 and AM 1146-270, are rather intriguing 
(even though we were only able to obtain an upper limit 
for AM 1146-270, it is still quite a restrictive limit). As dis- 
cussed earlier, the X-ray luminosities of these systems are 
low compared to ellipticals in general, more consistent with 
field spirals than with ellipticals. Furthermore, when com- 
pared to ellipticals of the same optical luminosity, it ap- 
pears that both galaxies may be relatively X-ray dim. This, 
combined with the fact that the ROSAT PSPC image of 
NGC 7252 (Fig. [n]) shows little in the way of any extended 
halo emission, poses a problem. 

In the previous two subsections, we have discussed the 
massive amounts of hot gas seen in the vicinity of these ear- 
lier interacting systems. The gas which is present in winds 
will leave the system on a short timescale, so the fact that 
these X-ray structures are not seen at the NGC 7252 epoch, 
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NGC 4038/9 Arp 220 NGC 2623 NGC 7252 

logLx (erg s" 1 ) 40.61 41.15 40.57 40.73 

T x (keV) 0.36 0.50 0.20 0.43 

n e (cm- 3 ) 0.0049 x 0.0047 x l/^/fj 0.0069 x l/^/rj 0.0105 x 1/^/fj 

A%as (M Q ) 1.5x10°^ 8.5 x lO 9 ^ 2.4 x lO 9 ^ 1.8 x lO 9 ^ 

E th (erg) 3.1 x lO 57 ^ 2.5 x lO 58 ^ 2.7 x lO 57 ^ 4.4 x lO 57 ^ 

t ?ool (Myr) 2410^ 5590^ 2290^ 2590^ 

A/ cooi ( M y r_1 ) °- 62 L52 im °- 70 

texp (Myr) 59 90 82 38 



Table 13. Values of physical parameters for the diffuse gas (using a hemispherical bubble model; see text) seen in the NGC 4038/9 
corona, the Arp 220 southern feature, the NGC 2623 eastern feature, and the NGC 7252 halo. r\ is the filling factor of the gas. 



a Gyr or so later, is not too surprising. The gaseous corona 
seen in the Antennae however, is different (coronal struc- 
tures such as that within the Antennae may exist within 
Arp 220 and NGC 2623 as well, though limitations on spa- 
tial resolution and on numbers of counts make it impos- 
sible to say). The Antennae's coronal gas, given the fact 
that it appears largely bound, with r\ close to unity (Read 
et al. 1995), is likely to have a long cooling time, perhaps 
1—2 Gyr, and hence such a gaseous halo should remain un- 
til the NGC 7252 epoch. This starburst generated gas may 
be the origin, as discussed earlier, of the hot gaseous halos 
known to exist around many large ellipticals (e.g. Fabbiano 
1989). For this hypothesis to work, one would expect hot 
gaseous halos around merger remnants. Why then, given 
that there is a good deal of evidence that hot gas should 
remain around the system until after obvious morphologi- 
cal evidence of a merger has vanished, does NGC 7252, the 
prototypical merger remnant candidate, show little sign of 
a hot halo? 

The problem is not in reconciling the observations of 
the earlier systems with those of NGC 7252. As can be 
seen from Table [l3[ the coronae seen in the Antennae and 
in NGC 7252 are very similar in terms of amount of gas. 
It is true though that, as the Antennae merger progresses 
through the ultraluminous phase, more energy and mass is 
likely to be injected into the halo, making it uncomfortably 
large when compared to the NGC 7252 X-ray observations. 
On the other hand cooling is taking place, returning the hot 
coronal gas to the disc, and over the 1.5 Gyr or so between 
the Antennae phase and the NGC 7252 phase, even at its 
present (relatively low) cooling rate, 1O 9 M0 of material will 
have returned to the central system. The real problem is 
that, in order for the merger hypothesis to work, and in or- 
der for these merger remnants to resemble elliptical galaxies 
in the X-ray, the halos need to be far larger (~ 50— 100 kpc), 
more massive (~ 10 9 - 10 10 M Q ; Fabbiano 1989; Forman et 
al. 1993) and hotter (T « 0.7 - 1.2 keV; Matsushita et al. 
1994; Rangarajan et al. 1995) than that observed around 
the prototype remnant, NGC 7252. 

The solution may be that we are looking for elliptical- 
like X-ray properties in merger remnants of far too young a 
dynamical age. For instance, it is thought that it will still 
take a few billion years for the optical luminosity, bright- 
ness distribution, colour gradients, velocity dispersion, UBV 
colours and fine structure (shells, loops, ripples etc.) of 
NGC 7252 to evolve to that consistent with a normal el- 



liptical (Schweizer 1996) . Perhaps it will take as long for the 
X-ray properties to do the same. 

In remants such as NGC 7252, quite a substantial 
amount of gas still remains, both in the tidal debris and in 
the central molecular disc. Hibbard & van Gorkom (1996) 
believe that NGC 7252 needs to get rid of half of its molec- 
ular gas and about 90% of its tidal material, in order to 
appear as an elliptical galaxy. The removal of the central 
molecular gas appears to present no real problem. What re- 
mains can be converted into stars within the next few Gyr, 
given current central star-formation rates (1— 2Mq yr~ ; Hi- 
bbard et al. 1994). The removal of the tidal debris however 
may be of particular interest. In a recent N-body simula- 
tion of NGC 7252 (Hibbard & Mihos 1995), it was seen that 
the majority of the material within the interaction-induced 
tidal tails is actually bound to the remnant. The bases of 
the tails fall back quickly, whereas the more distant ma- 
terial falls back more slowly. Some fraction of the material 
however may remain far from the remnant in very long tails, 
for very long times (e.g. 20% of the tidal material is believed 
to remain beyond 90 kpc for over a Hubble time). Although 
this material would be far too faint to observe, it may con- 
tribute to the Lya forest and to the IGM (Hibbard & van 
Gorkom 1996) . This return of the tidal debris therefore, can 
stretch out the merger remnant lifetime to many Gyr. 

The fact that the atomic gas disappears rather abruptly 
at the base of NGC 7252's northwestern tail, a region where 
the H I gas is seen to be falling toward smaller radii, whereas 
the optical tail continues across the face of the galaxy, is 
very suggestive of a quick and efficient process taking place 
converting the infalling gas to other phases (Hibbard et al. 
1994). One possibility as to what happens to the Hi gas is 
that it is converted into molecular gas through compression, 
in which case the system will remain rich in molecular gas 
for a long time. The other possibility is far more attractive; 
this returning material may be shock heated by the hot, 
though at this epoch, small halo. As this process continues, 
the X-ray luminosity and size of this halo will grow. The 
accompanying injection of energy will also raise the temper- 
ature. This is perhaps how large, hot X-ray halos are formed 
around merger remnants, (a requirement for the success of 
the merger hypothesis). The main bulk of the halo is not 
formed through the ejection of hot gas from the merger dur- 
ing the starburst stages of the interaction, but is instead 
formed, much later, through the shock-heating of tidal tail 
material falling back onto the remnant. If this picture is cor- 
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rect, then it is not surprising that only a small X-ray halo 
exists around NGC 7252. There is still a long time to wait (a 
few Gyr) before the bulk of the tail material returns, forming 
the large X-ray halos seen around many ellipticals. 

5.5 X-rays from tidal tails 

There is much evidence for enhanced star-formation, and 
even dwarf galaxy formation, within merger induced tidal 
tails (e.g. Schweizer 1978; Schombert et al. 1990). It is 
known that large amounts of neutral hydrogen are found 
within such tails (e.g. van der Hulst 1979; Hibbard et al. 
1994), and the self-gravity of this material can cause clump- 
ing, and lead eventually to star-formation. At the tip of the 
Antennae's southern tail, for instance, Mirabel, Dottori & 
Lutz (1992) have discovered a dwarf galaxy in the process of 
forming. This feature however, a chain of nebulae ionized by 
very young (< 6Myr) stars, is found to emit no significant 
X-ray emission (Read et al. 1995). This is not too surpris- 
ing, given that, if all the stars have been formed within the 
last 6 Myr, even the most massive of these would not have 
had time to evolve off the main sequence into the X-ray 
bright supernova remnant or high-mass X-ray binary phase. 
In later systems, however, X-ray features may be visible, as 
stars may have had time to evolve to X-ray bright phases. 

In practice, no significant emission is seen coincident 
with any tidal tail features anywhere within our merging 
galaxy sample. The soft diffuse feature seen to the west 
of NGC 2623, at first glance, may appear to have some- 
thing to do with the western tail. However, as discussed 
in Section ^(j, this is very unlikely to be the case. Even 
in NGC 7252, a system showing several large clumps of 
gas within its tidal tails (Hibbard et al. 1994), no obvi- 
ous X-ray counterparts are observed. A small amount of 
emission is seen however, within the eastern tail, spatially 
coincident with a giant Hn region (Hibbard et al. 1994) - 
the bluest region within the whole of NGC 7252. This fea- 
ture (the 'plume' of emission to the north of the source at 
a = 22ft20m56.25s, 8 = -24d42m06.4s) contains only ~ 7 
counts however, and a 2a upper limit to its (0.1— 2.0keV) 
luminosity (assuming a 5keV bremsstrahlung model) is 
Lx < 1.7 x 10 39 erg s _1 . Converting the optical luminosity 
of this Hn region (1.8xl0 41 erg s" 1 ; Hibbard et al. 1994) 
into an X-ray luminosity, assuming a similar Lx/Lb ratio 
to NGC 7252 as a whole, gives Lx = 1.7 X 10 38 erg s~\ a 
factor of ten below our upper limit. It is not, therefore, sur- 
prising that we do not detect any significant X-ray emission 
from this region. 



6 CONCLUSIONS 

We have completed a programme of ROSAT PSPC and HRI 
observations of a carefully chosen chronological sequence of 
eight merging systems in order to study the evolution of 
their X-ray properties through the merging process. While 
we have discussed the X-ray emission from each individual 
system in some detail, our main conclusions pertain to the 
evolutionary sample as a whole, and these conclusions are 
outlined below: 

• Across the sequence, from separated discs to relaxed 



remnants, the activity level is seen to rise, from a level be- 
tween that of nearby normal and starburst spirals, to a level 
far in excess of that of nearby starbursts, returning, at the 
end of the sequence to a level similar to that at the start. 

• Large amounts of hot gas are seen being forcibly ejected 
from these systems, the quantities of gas, and their extents, 
increasing with activity. The gas ejection process starts soon 
after the galaxies first meet each other, and evolves into what 
appear to be bipolar 'blowout'-type winds. These 'bipolar' 
features may be similar to those seen in nearby starbursts, 
the distortions due to the rapid evolution of the systems, 
or they may be due to two separate 'unipolar' winds, one 
ejected from each of the two interacting galaxies. What ap- 
pear to be massive, unipolar, ejections of hot gas are ob- 
served at the ultraluminous phases of the interaction, at the 
point of nuclear coalescence. These may result from a far 
more violent starburst taking place at this phase than at 
earlier stages, combined with very fast and vigourous dy- 
namical evolution of the system. At later times however, in 
the elliptical-like, relaxed remnant phase, almost all of this 
hot gas appears to have vanished, leaving only the sugges- 
tion of a small X-ray emitting halo. 

• The X-ray luminosity evolution, like the general activ- 
ity level, rises and falls along the evolutionary sequence. 
However, the amplitude of this change is considerably (more 
than an order of magnitude) smaller than that seen in the 
far-infrared. The resulting large drop in the Lx/Lfir ra- 
tio is probably realted to the ejected hot gas. The FIR flux 
is primarily due to dust heated by massive stars within the 
starburst, and hence tracks the star formation history. How- 
ever, the fraction of this energy which is radiated by hot gas 
is reduced as the starburst proceeds, as much of the input 
energy from supernovae and stellar winds actually goes into 
the kinetic energy of the ejected gas. 

• The prototypical merger remnant NGC 7252, bears lit- 
tle resemblance to typical elliptical galaxies in the X-ray, 
a fact at first appearing to conflict with the 'merger hy- 
pothesis', whereby elliptical galaxies can be formed from the 
merger of two disc galaxies. One possible resolution of this 
problem is that we are looking at far too early an evolution- 
ary age. Although NGC 7252's X-ray halo is far smaller, 
less massive, less bright and cooler than those of typical el- 
lipticals, shock heating of the returning tidal debris over the 
next few Gyr, may increasing the size, mass, luminosity and 
temperature of the halo. This may be the main source of an 
eventual elliptical-like X-ray halo, rather than the remnants 
of hot gas ejected into the halo at earlier phases. 

• The observations of NGC 520 reveal that it has no dis- 
cernable diffuse component, is very underluminous in the 
X-ray, and sits very uncomfortably in the X-ray evolution- 
ary sequence. It appears from other multiwavelength studies, 
that the NGC 520 merger has left the disc intact, whether 
because of gas-poor progenitors or a low progenitor mass 
ratio, and is likely not to evolve into an elliptical galaxy. 
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